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CHAPTER  I 


INTRODUCTION 

1 . 1  Introduction  and  Previous  Investigation 
The  sound  generated  by  turbomachines  is  not  only  of  high  ampli¬ 
tude  relative  to  the  ambient  background  noise  levels,  but  generally  its 
frequency  content  is  in  the  range  to  which  the  ear  is  particularly 
sensitive.  At  the  present  time,  the  increasing  site  and  wider 
application  of  such  machines  has  led  to  increased  interest  in  the 
reduction  of  turbomachine  noise. 

For  example,  the  noise  generated  by  a  jet  engine  involves  various 
acoustic  components:  externally,  the  exhaust  gas  mixing  process,  and 
internally,  the  unsteady  flow  through  the  fan  stage,  compressor,  and 
turbine.  A  significant  part  of  the  engine  noise  is  the  discrete 
frequencies  generated  by  the  interaction  of  the  rotating  and  stationary 
fan  blade  rows  at  the  funaamental  and  harmonics  of  the  fan  blade 
passing  frequency  (BPF).  Currently  the  BPF  generated  by  the  fan  in 
high-bypass  ratio  engines  dominates  the  acoustic  signature  of  the  jet 
engine.  However,  as  the  noise  generated  by  the  fan  is  reduced  with 
acoustic  treatment  in  the  fan  duct,  the  tone  noise  produced  by  the 
power  turbine  becomes  increasingly  important.  Because  the  mechanism 
of  BPF  generation  is  the  same  for  both  engine  components  the  study  of 
the  noise  generating  mechanisms  of  the  fan  will  lead  to  a  quieter 
engine.  This  is  one  reason  why  the  study  of  fan  noise  is  of  interest. 


The  emission  of  fan  noise  can  be  divided  into  three  phases: 

1 .  Generation 

Fan  noise  can  be  characterized  in  terms  of  two  component 
broadband  noise  of  relatively  low  amplitude  and  hip, her  ampli¬ 
tude  discrete  blade  passing  frequency  (BPF)  tones.  Broadband 
noise  results  mainly  from  (i)  the  random  fluctuation  in 
pressures  on  the  rotor  and  stator  components,  and  (ii)  the 
turbulent  inlet  flow  to  the  fan.  The  mechanism  that  generate 
the  discrete  BPF  tones  is  the  periodic  interaction  of 
downstream  fan  components  with  velocity  deficiencies  visoousl 
generated  hv  upstream  fan  components.  The  rotor  and  stator 
components  are  arranged  in  the  form  of  a  row  of  airfoils  know 
as  a  cascade.  F.ach  blade  in  the  cascade  is,  in  general, 
subject  to  disturbances  from  the  upstream  blades  as  well  as 
the  blades  in  the  same  row.  The  unsteady  lift  force  generate 
on  each  blade  is  equivalent  to  an  acoustic  dipole  when  the 
acoustic  wavelength  of  tho  unsteady  lift  is  much  gre;.'  or  than 
the  blade  chord.  When  this  condition  prevails,  tho  unsteadv 
lift  corresponds  to  a  "compact"  source.  A  compact  acoustic 
dipole  is  analogous  to  a  pair  of  pulsating  spheres  (point 
sources)  separated  by  a  small  distance  in  comparison  to  the 
acoustic  wavelength.  The  pulsation  of  the  spheres  is  in 
ant  i - phase  and  the  eonhinat  ion  generates  a  force  on  the  fluid 
along  a  line  connecting,  the  center  of  the  spheres.  I  i  1  ;  , 


f 1  no tunt ing  force  corresponds  to  the  unsteadv  lift  vector. 
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2 .  Transmission 

The  transmission  of  an  acoustic  wave  in  a  duct  can  be 
described  in  terms  of  acoustic  modes  (m,n).  A  specific  mode 
is  one  of  an  infinite  series  of  solutions  of  the  wave  equation 
for  the  specified  duct  geometry  and  boundary  conditions.  If 
the  acoustic  pressures  at  each  blade  in  a  cascade  are  in  phase 
with  each  other,  the  acoustic  energy  will  be  transmitted  along 
the  duct  in  the  plane  wave  (0,0)  mode.  This  node  has  a  uni¬ 
form  distribution  of  acoustic  pressure  over  the  duct  radius 
and  all  frequencies  propagate  axially  along  the  duct  at  the 
f roe-field  sound  speed.  It  is  also  possible  to  transmit 
acoustic  energy  along  a  duct  in  higher  order  modes.  These 
higher  order  nodes  have  a  nonuniform  spatial  distribution  of 
acoustic  pressure  and  above  the  respective  cutoff  frequency 
have  a  phase  speed  greater  than  the  free-field  sound  speed. 
Higher  order  modes  consist  of  non-spinning  modes  (0,n)  and 
spinning  modes  (m,n),  where  m  is  the  number  of  pressure  nodes 
in  the  radial  direction.  The  wave  front  of  non-spinning  modes 
lias  a  purely  axial  propagation  direction,  but  the  wave  front 
of  spinning  modes  (m,n)  travels  in  a  helix  along  the  duct.  Th 
solution  of  the  wave  equation  indicates  that  these  higher  mode 
will  decay  exponentially  along  the  duct  in  the  axial  direetion 
'.unless  the  exciting  f requetnev  is  greater  than  the  specific 
modal  cutoff  frequency.  Figure  1  shows  the  eh. tract  er  i -.t  ics  of 
the  transmission  of  the  plane  and  spinning,  acoustic  duet  mode; 

Tn  addition  to  the  plane  wave  mode  this  thesis  is 


specif icall'  concerned  with  the  first  higher  order  acoustic 
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mode  (-1,0).  The  minus  indicates  that  the  mode  is  spinning 
in  a  circumferential  direction  opposite  to  the  rotational 
direction  of  the  rotor.  The  one  indicates  there  is  one 
sinusoidal  cycle  of  pressure  variation  circumferentially 
around  the  duct.  The  zero  indicates  that  there  are  no  rever¬ 
sals  of  the  mode  pressure  along  a  radius.  A  further  property 
of  all  higher  order  modes  is  a  maximum  pressure  on  the  outer 
wall . 

3 .  Radiation 

For  tiie  propagating  modes,  there  is  a  strong  pressure 
field  at  the  open  end  of  the  duct  and  this  acoustic  pressure 
radiates  into  free  space  in  a  directivity  pattern  dependent 
upon  the  particular  mode.  For  the  plane  wave  mode,  the  far- 
field  acoustic  pressure  is  a  maximum  on  the  duct  axis.  In 
contrast,  for  the  propagating  (1,0)  mode  at  cutoff  frequency 
ratios  greater  than  one,  cancellation  occurs  on  the  duct 
centerline.  For  an  acoustic  dipole,  the  pressure  from  the 
positive  and  negative  sources  cancel  on  the  plane  normal 
to  the  dipole  axis.  The  directivity  pattern  on  the  (1,0) 
mode  at  frequencies  slightly  above  cutoff  resembles  (in  a 
plane  passing,  through  the  duct  centerline)  that  of  an  acoustic 
dipoli'  (lying,  in  this  plane)  witli  its  axis  oriented  across  the 
duct  opening. 

L.2  Previous  Theoretical  Work 

"lie  theorot  ical  ana  Ivr;  is  of  the  noise  from  turbomuch  i  lies  is  con- 
cern-s!  with  the  prediction  of  the  unsteadv  aerodvnanic  lift  and  tie- 
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far-field  pres  .tire  as  a  function  of  various  desiui  parameter.. 
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1.2.1  Unsteady  Lift 

The  interaction  of  a  cascade  of  airfoils  with  a  spatiallv  varying 
inlet  flow  is  shown  in  Figure  2.  The  velocity  deficiencies  of  the 
ups t roan  flow  cause  the  magnitude  as  well  as  the  direction  of  the 
velocity  relative  to  the  airfoils  to  change  with  tine,  thus  generating 
an  unsteady  lift  on  the  airfoils.  The  response  of  an  unsteady  flow  was 
first  solved  by  Ke:np  and  Sears  [!)*,  who  considered  an  isolated,  thin, 
flat-plate  airfoil  at  zero  incidence  in  two-dimensional  incompressible 
flow  subject  to  snail  velocity  nerturbat  ions  that  were  normal  to  the 
airfoil  chord.  later,  Horloek  12)  analvzed  the  effect  of  the  distur¬ 
bances  both  normal  and  parallel  to  the  airfoil  chord  on  an  isolated 
flat-plate  airfoil.  In  a  cascade  of  airfoils,  the  interference  between 
the  blades  in  a  row  becomes  espcciai.lv  s  i  pit  i  f  leant  for  high  solidity 
cases.  Whitehead  [3,  4]  analyzed  this  effect  for  a  single  infinite 
two-d inens ional  cascade  in  inviseiu  i ncoripross ibl e  flow.  An  analysis 
by  Henderson  and  Duneshyar  [  r> )  includes  the  effect's  of  blade  camber 
and  angle  of  incidence  in  the  cascade  response. 

Other  work  lias  been  performed  which  includes  the  compress i b i I  i  t v 
of  the  fluid  so  that  the  generated  sound  waves  can  he  calculated.  As 
an  example.  Smith  [ n j  considered  a  cascade  of  flat  plates  at  zero 
mean  incidence  to  he  replaced  hv  a  series  of  continuous  sin./.ul.iritv 
distributions  and  obtained  a  closed-form  solut  ion  lor  the  nnsfoadv 
lift  and  acoustic  pressure  due  to  several  perturb. i!  ion  tvs  ..  Another 
ex.'ifs  1 1  is  tiie  work  o:  Ka  j  i  and  I'k.izak  i  I  ’],  who  d  pressnt  e  d  t 

to  r  ■present  a  fiat  -plate  cascade  and  treated  it  hv  unit  led  1  i  r  <a  t  i  .•  •  . : 

in  Im.ic’m  !■  b  ;i-paa!e  He ‘  . a  cnees 


J 


equations.  Mani  [8]  assumed  the  blade  chord  and  spacing  to  be  small 
compared  with  the  wavelength  and  thus  simplified  the  cascade  response 
problem  to  that  of  an  actuator  disk  model.  Osborne  [9]  used  the 
linearized  potential  flow  equation  to  predict  the  unsteady  lift  and 
moment . 

1.2.?  Acoustic  Generation 

Tyler  and  Sofrin  [10]  studied  the  acoustic  excitation  and 
transmission  of  sound  from  a  known  pressure  source  in  an  annular  duct. 
They  demonstrated  that  aerodynamic  interaction  between  blade  rows 
could  result  in  a  rotating  pressure  pattern.  This  rotating  pattern 
could  excite  propagating  spinning  modes  which  require  a  supersonic 
tangential  phase  speed,  even  though  the  rotor  tip  velocity  is  subsonic 

1.2.3  Acoustic  Transmission 

The  description  of  the  transmission  of  sound  in  an  infinite  duct 
in  terms  of  Bessel  and  Neumann  functions  has  been  known  lor  some  time. 
More  recent  work  lias  been  performed  by  Morfev  [11],  who  investigated 
the  reflection  and  radiation  of  sound  from  the  end  of  a  duct.  The 
transmission  of  sound  from  the  end  of  the  duct  is  the  radiation  probie 
and  depends  on  the  boundary  conditions  at  the  duct  exit.  For  example, 
Tyler  and  Sofrin  [10]  used  a  baffled  duct  model  to  predict  tin  aconst i 
measurements  of  a  fan  for  angles  not  near  90  degrees. 

1.3  1’revious  l'.xpe  r  i  non!  a  1  Work 

per  imer.l  a  1  work  has  concent  rat  ed  on  measuring,  tin  unsteady  liit 
res;-  ir.se  of  a  blade  in  a  cascade  to  measured  inflow  distortions  withou 
me  i.-i  rotten  t  ot  the  acoustic  output.  In  those  instances  where  acoustic 


t 


Q 


pressures  have  been  measured,  only  the  inflow  distortion  lias  heen 
measured;  the  unsteady  lift  has  not  been  sinml taneousl v  measured.  The 
uns toady  lift  response  of  a  blade  in  a  cascade  has  been  measured  by 
Bruce  [12].  Recently,  measurement  of  the  fluctuating  pressure  on 
blades  has  become  possible  with  the  development  of  miniature  transducers 
Satyanarayana ,  Henderson,  and  Gostelow  [13]  used  'l’itran'  pressure 
transducers  inserted  in  an  uncambered  stator  of  6-inch  chord  in  a  single 
stage  fan.  They  measured  the  pressure  on  both  sides  of  the  blade  at 
15  points  along  the  chord.  These  measurements  were  found  to  be  in 
agreement  with  the  analysis  of  Reference  [4].  Callus  et  al.  [1 4J  used 
'Kulite'  transducers  to  measure  the  fluctuating  pressures  on  both 
rotor  and  stator  blades  along  the  chord  at  the  midspan. 

Acoustic  measurements  of  the  acoustic  tonal  energy  generated  by 
a  fan  have  been  performed  be  several  investigators.  Early  work  in  this 
area  is  that  of  Shari  and  [15]  who  found  that  the  discrete  frequency 
noise  from  multi-stage  fans  arose  from  the  aerodynamic  interaction  of 
the  fixed  and  rotating,  blade  rows.  He  also  reported  that  the  amplitude 
of  this  tone  noise  was  strongly  dependent'  upon  the  axial  spacing  between 
the  blade  rows. 

Other  early  work  is  that  of  Bragg  and  Bri  Ige  [16]  who  report  the 
re. nils  of  a  panel  of  experts.  They  concluded  that  the  Lone  noise 
from  turbojet  compressors  was  dipole  in  origin  and  predominately  from 
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two  BPF  harmonics,  and  obtained  a  reasonable  agreement  with  the  theory 
of  Reference  [8].  Hard  and  Perulli  {20]  employed  a  six-bladed  fan 
at  low  axial  flow  speed  and  measured  the  frequency  spectra,  axial , 
circumferential,  and  radial  wave  number  of  the  sound  field  in  order 
to  compare  the  structure  of  the  acoustic  modes  in  the  duct  with  theory. 
Snith  [6]  measured  the  inflow  distortion  generated  by  screens  and 
resulting  acoustic  pressures  in  a  high  hub-to-tip  radius  ratio  rotor. 
After  accounting  for  reflection  at  the  duct  exit,  he  obtained  excellent 
agreement  with  his  theory.  Moore  [21]  used  a  microphone  cross-correlated 
with  a  reference  signal  to  measure  the  rotating  sound  field.  Yardley 
[22]  proposed  a  wave  number  spectrum  measurement  technique  to  analyze 
the  acoustic  field  in  a  duct. 

The  possibility  that  an  acoustic  field  can  strongly  influence  an 
essentially  incompressible  hydrodynamic  field  was  suggested  by  Doak  [23]. 
He  called  this  effect  an  acoustic  back-reaction.  Archibald  [24,  25]  in 
two  companion  papers  found  that  indeed  a  high  amplitude  acoustic  field 
was  able  to  control  the  two  instability  phenomena  of  vortex  shedding 
and  Tol lm ion-Schl i cht ing  waves.  Crow  and  Champagne  [26]  found  that  an 
acoustic  field  was  able  to  produce  a  coherent  structure  in  the  flow 
from  a  low  Mach  number  jet. 

In  addition  to  those  references  which  deal  directly  with  the 
problem  order  investigation,  Morphey  [27]  and  Cumpstey  [ 2S ]  have  written 
critical  reviews  of  engine-related  aero-neons t ic  noise  souces. 

1.4  Present  Work 

!n  .u  rv-'v  i  ny  tie-  existing,  knowledge  of  this  field,  it  appears  that 
there  ate  sev-Til  prediet  ion  formulae  for  unsteadv  lift  as  well  as 


acoustic  modal  pressures,  and  the  experimental  work  is  in  good  agree¬ 
ment  with  the  theory.  But,  simultaneous  detailed  measurements  of  the 
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aerodynamic  pressure  on  the  blades  and  the  acoustic  modal  pressure  in 
a  duct  are  necessary  to  confirm  the  theoretical  prediction  and  to 
find  the  coupling  coefficient  of  the  unsteady  lift  to  acoustic  modal 
pressure.  In  order  to  remove  the  possibility  of  ambiguity  in  this 
matter,  the  present  work  has  been  undertaken  to  generate  and  measure 
the  acoustic  pressure  of  duct  modes  due  to  the  interaction  mechanism, 
and  simultaneously  to  measure  the  unsteady  lift  on  a  rotor  blade.  This 
will  establish  the  coupling  factor  which  has  not  received  attention 
previously,  and  will  permit  investigation  of  any  back-reaction  effects 
of  the  acoustic  pressures  on  the  unsteady  aerodynamic  response. 


CHAPTER  II 


EXPERIMENTAL  ASPECTS 

2. 1  Introduction 

The  experiment  was  designed  to  investigate  the  generation  of  fan 
blade  passing  frequency  (BPF)  tone  noise  in  a  duct.  Figure  3  shows 
the  theoretical  mechanisms  and  the  measurements  that  were  performed  at 
various  positions  in  the  noise  generation  process.  These  measurements 
were  used  to  confirm  the  theory  at  many  points  in  the  generation 
process.  The  velocity  field  was  produced  by  distortion  screens  and 
measured  using  a  five-hole  probe.  The  circumferential  velocity  profile 
was  then  Fourier  analyzed.  The  rotor  blades  interacting  with  the 
distorted  velocity  profile  generate  an  unsteady  lift  which  excites  acous¬ 
tic  modes  in  the  duct  and  eventually  radiate  as  BPF  sound.  The  unsteady 
force  normal  to  the  chord  (lift)  was  measured  using  a  strain  gage 
sensor  installed  inside  a  one-inch  spanwise  segment  at  midspan  of  a 
rotor  blade  shown  in  Figure  4.  The  acoustic  pressure  was  measured 
with  a  flush-mounted  1/8-inch  Bruel  and  Kjaer  (B&K)  microphone  at  four 
locations  along  the  duct  to  decompose  the  pressure  into  duct  modes  and 
to  find  the  coupling  coefficient  between  the  aerodynamic  pressure  and 
the  acoustic  modal  pressure. 


2 . 2  Apparatus 

The  experiment  was  performed  in  the  Axial  Flow  Research  Fan  (AFEF) , 
a  facility  of  the  Garfield  Thomas  Water  Tunnel  which  is  located  at  The 
Pennsylvania  State  University.  This  facility  is  shown  in  Figure  3. 
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This  apparatus  has  been  described  by  bruce  [29].  The  AFRF  is  19.6 
feet  long  and  consists  of  a  bellmouth  inlet  leading  to  an  annular  flow 
passage,  which  contains  a  test  rotor  and  an  auxiliary  fan.  At  the  end 
of  the  AFRF,  there  is  an  exhaust  throttle.  The  facility  has  a  hub 
radius  of  A. 75  inches  and  a  tip  radius  of  10.75  inches.  The  hub-tip 
radius  ratio  is  0.442. 

The  auxiliary  fan  is  driven  by  a  70  HP  motor  and  is  capable  of 
delivering  1,500  cubic  feet  of  air  per  minute  at  a  pressure  of  3.5 
inches  of  water  gage.  The  research  rotor  is  driven  by  a  separate  20  HP 
motor;  thus,  in  conjunction  with  the  auxiliary  fan,  the  research  rotor 
can  be  operated  over  a  range  of  steady  lift  coefficients.  The  research 
rotor  consists  of  nine  cambered  blades  with  a  chord  length  of  six  inches 
and  a  span  of  six  inches.  The  blade  profile  was  a  circular-arc  camber- 
line  with  10  percent  maximum  thickness-to-chord  ratio.  A  list  of 
rotor  parameters  is  given  in  Table  1  and  the  details  of  the  rotor 
design  are  described  in  Appendix  A. 

The  four  microphone  measurement  positions  were  selected  to  optimise 
the  discrimination  between  the  (-1,0)  and  (0,0)  modes  after  considering, 
acoustic  reflection  effects  from  the  duct  open  end  and  t he  rotor.  The 
details  of  the  selection  of  the  microphone  locations  Is  described  in 
Appendix  b.  Figure  6  shows  the  locations  of  the  distortion  screen, 
test  rotor,  microphone;;,  and  static  pressure  tap  at  the  inlet  for 
dot  urn i n i ng  the  mean  axial  velocity  V  .  The  design  procedure  of  the 
distortion  screens  is  presented  in  Appendix  ll. 
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2 . 3  last  rumen t at i on 

A  five-holi-  probe  located  in  the  plane  of  the  rotor  leading  edge 
was  employed  to  measure  the  local  velocity  with  the  rotor  removed.  The 
pressures  from  each  hole  in  the  probe  were  converted  into  DC  voltages 
by  using  five  separate  pressure  transducers.  The  voltages  were  proces¬ 
sed  by  a  multi-channel  integrating  digital  voltmeter  with  the  output 
punched  onto  a  paper  tape.  The  paper  tape  was  then  fed  into  a  computer 
for  data  reduction.  The  measured  velocity  profiles  were  Fourier 
analyzed.  The  magnitude  and  phase  of  30  harmonics  were  determined.  A 
block  diagram  of  the  instrumentation  is  shown  in  Figure  7. 

A  Harrison  amplifier  was  used  to  amplify  the  voltage  of  the  AC 
signal  from  the  instrumented  rotor  blade  by  1,000  times  to  supply  a 
constant  tour  volts  DC  to  the  strain  gage  through  the  slip  rings 
inside  the  hub.  A  photodiode  sensed  the  light  passing  through  a  single 
slit  in  a  disk  on  the  rotor  shaft  generating  a  onee-per-revolut ion 
pulse  as  a  c.  i  rcumf  vrent  ial  position  reference.  Another  photodiode 
sensed  the  light  passing  through  an  inner  row  of  60  slots  on  the  same 
disk.  This  gave  a  direct  reading  of  the  shaft  speed  in  rpn  in  a  one- 
second  counting  period.  The  unsteady  pressure  was  measured  bv  one 
1/tf-incli  Hf.K  microphone  connected  to  a  li.Sk  2603  amplifier  which  also 
supplied  the  DC  voltage  to  the  microphone.  This  microphone  was  then 
moved  system. »t  i  cal  1  v  to  tile  other  locations  in  the  duct  which  were 
eleid  with  flush  tilting  inserts. 

iiie  si-ti.il'  iron  the  lift  gage  and  microphone  were  anal  y.-.e.i  using, 
a  ussier  «>:  the  .malvsi:.  options  of  the  Spectral  Dvnimies  Corporation 
1  hit)  t  wo  -  e  is  sine  1  digital  signal  processor.  Tin  SD-5M)  has  a 


ha;  It -in  f'a.t  Fourier  Tran  stot — >  (ITT)  which  can  produce  an  ensemble 
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averaged  spectrum  (function  A).  In  conjunction  with  a  Spectral 
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Dynamics  5D-134A  signal  racio  tuner,  it  can  perrons  phase-locked 
averaging  in  the  time  domain  (function  12)  and  phase-locked  averaging/ 

FFT  analysis  (function  13).  The  time  position  reference  is  the  once- 
per-revolution  pulse.  Figure  8  shows  a  schematic  diagram  of  this 
instrumentation.  The  tracking  ratio  tuner,  using  the  once-per-revolut ion 
pulse,  generates  exactly  1024  pulses  per  revolution  of  the  rotor.  These 
pulses  trigger  the  SD-360  to  sample  the  input  signal  and  completely  fill 
the  memory  of  the  analyzer  in  one  shaft  revolution.  Tine  samples  of  the 
signal  of  each  subsequent  shaft  revolution  are  added  to  the  average  of 
the  previous  data  in  the  memory  to  obtain  an  ensemble  average  in  the 
time  domain. 

The  raw  signals  from  the  strain  gage  and  microphone  were  also 
connected  to  a  two  channel  oscilloscope  for  a  visual  check  of  the 
quality  of  the  signals.  A  pressure  transducer  and  its  associated 
signal  conditioning  equipment  ver?  used  to  measure  the  wall  static 
pressure  of  the  duet  for  calculating  the  mean  axial  velocity  for 
each  test. 

The  setting  of  the  frequency  rnn.e  and  aliasing  filter  on  the 
signal  processor,  and  the  frequency  ratio  on  the  trackin';  ratio  tuner 
are  important  for  the  time  domain  averaging  ( functions  12  .rad  13).  A 
description  oi  the  effect  of  variation  in  these  settings  is  contained 
in  Appendix  D. 
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2 . 4  Calibrnt i  on 

2.4.1  Five-hole  Probe  Calibration 

The  calibration  of  the  probe  involved  experimentally  determining 
a  set  of  data  in  a  known  flow  field  that  was  processed  by  a  computer 
and  resulted  in  a  description  of  the  response  of  the  probe.  The  probe 
was  mounted  in  a  calibration  device  which  permitted  a  +30-degree 
rotation  in  both  the  yaw  and  pitch  planes.  The  pressure  coefficient  as 
a  function  of  flow  velocity  and  angle  was  measured  in  the  potential 
core  of  an  open-jet  facility  at  each  calibration  point.  A  detailed 
description  of  the  calibration  procedure  was  presented  by  Troaster  and 
Yocum  [30].  The  calibration  curve  for  the  probe  used  in  this  investi¬ 
gation  is  shown  in  Figure  9. 

2.4.2  Strain  Gage  Calibration 

A  static  calibration  of  the  strain  gage  was  performed  by  varying 
the  masses  loaded  on  the  instrumented  blade;  this  resulted  in  a  sensi¬ 
tivity  to  a  normal  force  of  0.972  lb/volt.  The  dynamic  response  of  the 
gage,  shown  in  Figure  10,  was  obtained  from  the  broad  band  envelope  of 
the  spectrum  of  the  gage  signal  when  no  distortion  screen  was  installed 
The  resonance  of  the  blade  section  was  found  to  be  on  the  order  of  3(10 
!!.:  which  is  well  above  the  cutoff  frequency  of  frequency  of  the  (1,01 
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microphone.  The  spectrum  with  the  piston  phone  applied  is  shown 
in  Figure  11  and  resulted  jn  a  57.64  dB  voltage  level  (reference  to 
0.1  volts  peak)  with  the  B&K  amplifier  and  SD-360  analyzer  attenuators 
at  zero  for  this  end-to-end  calibration. 


2.4.4  Signal  Processor  Calibration 

The  SD-360  provides  an  internal  calibration  tone  of  0  dB  or  676 
linear  scale  units,  relative  to  0.0707  volts  rms  at  approximately  25 
percent  of  the  full  frequency  range.  The  0  dB  reference  was  set  using 
this  calibration  tone  every  time  the  power  switch  was  turned  on  or  a 
different  signal  analysis  option  was  selected.  A  chock  of  the  internal 
calibration  was  also  conducted  using  a  signal  generator  and  voltmeter 
and  it  was  found  that  0.071  volts  rms  input  sinusoidal  signal  corres¬ 
ponded  to  0  dB,  or  6S0  linear  scale  units.  Tills  is  the  same  as  the 
internal  calibration  within  the  experimental  accuracy. 

The  calibration  of  the  strain  gage  and  the  signal  processor  was 
also  checked  by  using  the  effect  of  gravity  on  the  instrumented  blade 
segment.  The  gravitational  force  produces  a  one-cvcle  sinusoidal  signal 
on  the  lift  gage  in  one  revolution  of  the  rotor  blade.  Therefore,  tonal 
energy  at  the  fundamental  (i.e.,  shaft  rate)  in  the  first  spectrum  of 
the  lift  gage  signal  should  he  the  result  of  the  gravitational  force  and 
be  independent  of  rotor  speed.  The  magnitude  of  this  harmonic  was  found 
to  be  0.033  and  0.060  volts  peak  for  all  rotor  speeds,  when  the  nine- 
and  Len-cvc 1 e  screens  were  installed,  respectively.  The  weight  of  the 
blade  segment  was  0.056  pounds.  Using  this  information,  a  calibration 
constant  of  0.966  or  0.933  lb/volt  was  obt  •lined  tor  the  two  separate 
ra.v.i.  This  is  a  good  confirmation  of  the  static  calibration  vain,  of 
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0.972  lb/volt.  Confirmation  of  the  validity  of  the  gravitational 
check  is  the  fact  that  the  phase  angle  of  the  first  harmonic  maintained 
a  constant  value  of  -73  degrees  for  all  rotor  speeds.  The  gravitational 
force  acts  downward  at  the  -90  degrees  position  from  the  12  o'clock 
location  of  the  duct  and  produces  a  postive  value  of  lift.  Since  the 
once-per-revolution  pulse  was  at  -16.2  degrees,  the  maximum  lift  should 
occur  at  (-90  +  16.2)  -  73.8  degrees.  Consequently,  the  effect  of 
gravity  is  apparently  producing  the  first  harmonic  on  the  lift  spectrum 
within  the  measurement  accuracy. 

2 . 5  Measurement 

The  disturbance  screens  were  located  at  16.5  inches  upstream  of 
the  leading  edge  of  the  rotor.  The  distance  between  screen  and  rotor 
is  determined  by  considering  the  effect  of  turbulent  mixing  downstream 
of  the  screens.  Turbulent  mixing  tends  to  reduce  large  velocity  grad¬ 
ients  caused  by  step  changes  in  the  flow  resistance  coefficient,  K,  of 
the  screens.  Far  downstream  of  the  screen  the  velocity  profile  becomes 
smoother,  however,  the  fundamental  harmonic,  Af.also  has  decayed,  which 
is  not  desired.  The  optimal  distance  is  larger  for  screens  with  low 
numbers  of  circumferential  distortion  cycles  than  screens  with  high 
numbers  of  distortion  cycles  because  the  velocity  gradient  is  smaller 
for  the  low  number  distortion  cycle  case  and  there  is  consequently  less 
turbulent  mixing. 

Measurements  were  conducted  to  find  the  characteristics  of  the 
flow  fields  produced  by  screens  as  a  function  of  circumferential  angle, 
3,  at  three  radial  positions:  midspan  (r  =  7.75  inches),  one  inch  from 
the  rotor  tip  (r  =  9.75  inches),  and  one  inch  from  the  rotor  .tub 
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(r  =  5.75  inches).  The  circumferential  surveys  were  performed  in 
2-degree  increments  by  keeping  the  screen  spatially  fixed  and 
rotating  a  downstream  portion  of  the  outer  casing.  The  casing  portion 
containing  the  velocity  probe  was  located  in  the  region  of  the  test 
rotor.  The  rotor  was  removed  during  the  velocity  surveys  to  avoid 
contaminating  the  measurement. 

The  mean  axial  velocity,  Vx,  was  selected  to  give  the  desired 

incidence  angle  over  the  operating  speed  range  of  the  rotor  of  the  AFRF. 

The  velocity  V’  can  only  be  calculated  when  the  velocity  profile  is 

known  and  it  is  not  practical  to  use  the  five-hole  probe  to  set  \’x. 

Thus,  tin'  static  pressure  AP  measured  with  a  wall  static  tap  just- 

downstream  of  Lhe  duc.t  inlet  when  the  screens  were  calibrated  was  used 

to  obtain  the  desired  velocity  in  the  following  manner.  Temperature 

and  atmospheric  pressure  were  recorded  to  calculate  the  air  density,  , . 

Using  Bernoulli's  equation,  the  \;x  can  he  related  to  the  pressure 

difference  AP  and  r.  The  quantity  AP s/p  plotted  versus  V  is  a  straight 

line  on  log-log  paper,  which  is  different  for  each  screen,  as  shown 

in  Figure  12.  Mien  V  =  0,  AP  / ;  represents  half  the  pressure  loss  due 

x  h 

to  tile  distortion  screen.  By  interpolation  of  these  plots,  the  ties  i  red 
value  of  \P  at  each  Vx  can  be  obtained  for  all  test  conditions.  The 
pressure  measured  from  the  Lap  was  then  used  for  set  Ling  the  operating 
condition  during  the  unsteady  pressure  measurements  portion  of  the  test. 

Tile  zero  e  i  reuml  isrent  i  a  1  position  was  selected  as  the  12  o'clock 
position  of  tin-  duet.  The  onee-por- revo 1 u t i on  pulse  was  measured  it' 
ore  i  r  .It  -lh. 2  degrees  relative  to  this  Zero  position,  i.e..  It'..'  dr'lere 

ot  rotor  rotation  prior  to  the  top.  Figure  1J  shows  the  spat  i.il  1  t  ion; 
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TABLE  II 


TEST  MATRIX  OF  UNSTEADY  PRESSURE  MEASUREMENTS 
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Figure  13  Spatial  Relation 
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used  in  the  phase  measurement  of  unsteady  lift.  The  phase  angle  of 
unsteady  lift  was  defined  relative  to  the  leading  edge  of  the  blade. 

When  the  blade  coincides  with  the  minimum  axial  velocity,  a  positive 
upwash  and  incidence  angle  is  defined.  The  unsteady  lift  signal  from 
the  lift  gage  was  electrically  positive  for  a  positive  deflection, 
which  is  in  the  same  direction  as  the  positive  upwash.  The  minimum 
axial  velocity  of  the  screen  was  positioned  at  the  12  o'clock  position. 
Note  that  one  cycle  of  the  nine-cycle  screen  only  occupies  l/9th  of 
the  360  degrees  in  the  duct,  i.e.,  40  degrees.  Thus,  the  16.2  rotor 
degrees  corresponds  to  (16.2/40  *  360)  145.8  degrees  in  the  screen  or 
unsteady  lift  cycle.  Similar  considerations  hold  for  the  ten-cycle 
screen. 

The  test  matrix  of  unsteady  pressure  measurements  is  shown  in 
Table  II.  The  variables  of  the  test  are  enumerated  below. 

1.  Number  of  cycles  of  the  disturbance  screen:  Nine-  and  ten- 
cycle  screens  were  used  to  generate  the  plane  wave  (0,0) 
mode  and  the  first  higher  order  (-1,0)  mode  for  the  nine- 
bladed  rotor.  Measurements  with  no  distortion  screen  were 
also  conducted  to  find  the  background  noise  level. 

2.  Research  rotor  speed:  Twelve  rotor  speeds  were  selected  to 
have  the  blade  passing  frequencies  (KPF)  below,  at,  and  above 
the  (1,0)  inode  cutoff  frequency  of  the  duct. 

3.  Flow  coefficient  <> :  The  mean  axial  velocity,  V'x,  was  set 

by  the  auxiliary  fan  so  that  the  flow  coefficient  6  =  V  /rui 

x 

could  be  selected  independently  of  the  research  rotor.  The 
mean  incidence  angle,  in,  at  midspan  was  determined  from 
velocity  triangles.  Two  flow  coefficients  $  "  0.70  ( i  -  0.3°) 
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and  <}>  =  0.56  (i  -  6.0°)  were  chosen  to  investigate  the 
effect  of  steady  load  on  the  unsteady  lift  and  acoustic 
pressure  measurements. 

2.6  Experimental  Accuracy  and  Repeatability 

2.6.1  Velocity  Survey 

In  the  velocity  measurements  downstream  of  the  distortion  screen, 
the  2-degree  circumferential  increment  of  the  outer  casing  was  obtained 
by  a  Slo-Syn  motor.  This  motor  was  controlled  by  an  indexer  capable 
of  a  minimum  of  1/80-degree  adjustment.  Thus,  very  little  error  would 
result  from  the  measurement  locations.  Measurements  were  repeated  four 
times  at  each  location;  a  statistical  computer  subroutine  was  then 
employed  to  analyze  the  data.  It  indicated  that  the  normalized  90 
percent  error  of  the  mean  for  the  voltage  reading  corresponding  to  the 
velocity  in  the  axial  direction  was  less  than  10  percent.  This  was  a 
typical  result  for  the  highly  turbulent  flow  field  behind  the  screens. 

The  amplitude  of  harmonic  numbers  N  between  20  and  30  were  found  to  be 
always  less  than  five  percent  of  the  fundamental  one. 

2.6.2  Unsteady  Lift  and  Pressure  Measurements 

The  raw  signals  from  the  strain  gage  and  the  microphone  were 
ensemble  avaraged  256  times  for  both  the  spectrum  and  the  phase-locked 
option.  The  phase-locked  average  (functions  12  and  13)  preserves  the 
signal  that  is  periodic  with  the  rotor  speed,  while  random  or  non-periodic 
signals  will  tend  to  diminish  due  to  the  large  number  of  sums  in  the 
average.  Thus,  the  magnitude  and  the  phase  obtained  were  due  to  the  peri¬ 
odic  interaction  between  the  spatially  fixed  velocity  distortion  and  the 
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rotor.  Furthermore ,  the  values  of  the  magnitude  obtained  with  the 
spectrum  and  phase-locked  average  were  in  close  agreement,  indicat  ini; 
that  there  was  no  significant  effect  of  inflow  turbulence  causing  high 
levels  of  intermittent  spatial  distortion. 

Each  unsteady  measurement  was  repeated  at  least  four  times  on  a 
different  date.  It  was  found  that  the  acoustic  pressure  repeated  very 
well;  the  spread  in  magnitude  and  phase  were  less  than  +0.2  dB  and 
+2  degrees.  The  differences  of  magnitude  and  phase  on  the  unsteady  lift 
were  within  +0.5  dB  and  +5  degrees  for  most  cases.  However,  larger 
fluctuations  of  lift  were  observed  for  speeds  near  the  cutoff  frequency 
about  +1-0  dB  and  +20  degrees  for  the  worst  case.  The  normalized 
standard  deviation  of  the  unsteady  measurement  was  generally  less  than 
4.5  percent.  It  increased  to  23.4  percent  for  the  largest  fluctuation 
of  lift  near  cutoff. 

2 . 7  Pat a  Validity 

Figure  14  is  a  spectral  comparison  of  the  acoustic  pressure  wit’a  the 
rotor  operating  above  cutoff  when  the  nine-  and  ten-cycle  screens  were 
installed.  The  background  noise  was  measured  at  the  same  shaft  speed, 
but  with  no  distortion  screen  installed.  Because  the  background  noise 
was  lower  by  30  dB  than  the  signal  at  BPF,  the  experimental  measurements 
are  considered  valid. 

Figure  15  is  a  spectral  comparison  of  the  unsteady  lift  gage  with 
the  rotor  operating  above  cutoff  when  the  ten-cycle  screen  was  installed. 
(The  ntne-cycle  screen  gave  similar  results.!  The  background  noise  was 
measured  at  the  same  shaft  speed,  but  with  no  distortion  screen  installed. 
Another  form  of  background  noise  is  the  electrical  noise  due  to  the  slip 


igure  15  Background  Uoise  of  Unsteady  Lift  Gage  Measurement 


rings.  This  was  determined  by  replacing  the  strain  gage  with  an 
equivalent  resistor  bridge.  The  output  of  this  configuration  at 
the  same  operating  conditions  is  also  shown  in  Figure  15.  Although 
the  tones  at  integer  multiples  of  the  shaft  rotation  frequency  were 
generally  the  same  with  or  without  the  screen  present,  a  25-dB  difference 
in  level  was  observed  at  the  tenth  harmonic  of  the  shaft  frequency. 

Thus,  the  signal-to-noise  ratio  is  sufficient  to  make  the  data  valid. 


CHAPTER  III 


RESULTS  AND  COMPARISON  WITH  THEORY 

3. 1  The  Sound  Field  in  a  Duct 

The  modal  structure  of  the  sound  field  in  an  infinite-length 
annular  duct  with  no  mean  axial  flow  has  been  well  studied  and  under¬ 
stood  [31,  32].  The  resulting  sound  pressures  can  be  represented  in 
the  form  of  an  array  of  transformed  modal  distribution  functions.  The 
circumferential  part  is  simply  a  Fourier  transformation  and  consists 
of  corotating  and  counterrotating  modes.  The  radial  transform  is 
more  complicated;  the  resulting  characteristic  functions  depend  on 
circumferential  mode  number,  hub-tip  radius  ratio,  and  the  radial 
mode  number.  The  transformed  wave  number  in  the  axial  direction 
determines  the  propagation  characteristics;  the  transition  from  acoustic 
pressure  decay  to  propagation  occurs  when  the  circumferential  Mach 
number  of  the  rotating  mode  reaches  a  critical  value.  Because  the 
axial  Mach  number  was  less  than  0.2,  the  effect  of  the  convection 
velocity  was  neglected  in  this  initial  investigation. 

The  sound  field  in  an  annular  duct  with  no  mean  flow  can  be 
expressed  as 

P(r ,0,x,t)  -  mEn  pnn  (kmnr)  exp  i  (a^  kmn  x  +  mO  -  mt),  (1) 

where  m  is  the  circumferential  mode  order,  which  can  be  positive 
(corotation)  or  negative  (counterrotation)  relative  to  the  direction 
of  rotation  of  the  rotor.  The  radial  mode  orcer  is  indicated  by  n. 


40 


The  coefficient  of  the  wave  number  in  the  axial  direction,  a  ,  is 

mil 

given  by: 


k 

x 

a  =  - - 

mn  k 

mn 


(2) 


where  y  is  the  cutoff  ratio,  defined  by  the  ratio  of  the  frequency 


of  the  rotating  mode  f  and  its  cutoff  frequency,  f*.  The  cutoff 


frequency  is  related  to  the  characteristic  wave  number  of  the  mode 


k  and  the  sound  speed  as: 
mn 

k  c 

f*  =  — 

2  7T 


In  equation  (1),  E  ^  is  the  radial  pressure  distribution  function 
formed  by  the  Bessel  function  of  the  first  kind,  J,  and  the  second 
kind,  Y;  i.e. : 


E  =  C  [J  (k  r)  +  Q  Y  (k  r)] 
mn  mn  m  mn  rain  m  mn 


(3) 


The  values  of  Q  and  k  are  determined  by  the  boundary  conditions, 
mn  mn 

The  radial  pressure  gradient  must  be  zero  at  the  inner  and  outer  walls, 
therefore : 


and 


J'(k>  )  + 
m  mn 


J'(k' 
m  mn 


o) 


Y 

mn 


(k’>  =  0 
mn 


Y'  (k'  o)  =  0 
ran  mn 


(4) 


where  k'  is  the  characteristic  wave  number,  k  ,  normalized  by  the 
ran  mn 

outer  radius,  r.  .  The  hub-tip  radius  ratio  is  o.  The  orthogonality 
b 
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r 

? 


of  the  characteristic  functions  requires  that 

r,  0  p,  v  m,  n 

b 

E  (r)  E  (r)  2tt  r  dr  =  if  .  (5) 

ran  pv 

r  1  p ,  v  =  m ,  n 

3 


The  square  of  the  modal  pressure  represents  the  power  of  the 

modes;  thus,  the  normalized  constant,  C  ,  ensures  that  the  integral 

mn 

2 

of  p  across  the  duct  is  unity.  Combining  Equations  (3)  and  (5) 

gives,  for  the  value  of  C  : 

mn 


2 

I —  “ 

2 

a  - 

m 

J  (k‘  o)  +  Q  Y  Ck ’  o) 

k* 

m  mn  mn  m  mn 

mn 

L  A 

_ 

_ 

(6) 


If  the  frequency  of  the  exciting  source  of  a  mode  exceeds  the 
cutoff  frequency  of  that  mode,  the  acoustic  pressure  propagates,  and 
the  wave  front  travels  at  a  helix  angle,  a,  relative  to  the  axial 
direction.  The  value  of  a  is  given  by: 

a  =  sin  ^  (1/y) 


For  the  plane  wave,  the  propagation  direction  is  purely  axial, 
i.e.,  a  =  0°.  For  higher  order  spinning  modes  just  at  cutoff,  the 
wave  fronts  rotate  in  the  circumferential  direction,  i.e.,  a  =  90°. 

The  value  of  the  characteristic  constants  for  the  AFRF  geometry 
are:  k'  Q  =  1.425;  Q1Q  =  0.258;  Q  =  1.018. 


% 
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3.2  Aerodynamic  Lift  bred  i  c  t  ion 
Whitehead  [3,  4]  represented  a  thin  zero-camber  cascade  by  a 
number  of  t nually  spaced  vortices  with  continuously  varying  strength 
to  obtain  a  solution  of  the  chordwise  pressure  distribution  for 
various  forms  of  perturbation.  Smith  [hi  modified  Whitehead's  model 
to  include  compressibility  effects.  The  unsteady  lift  results 
referred  to  later  in  this  thesis  are  obtained  bv  integration  of  the 
unsteady  pressures  obtained  from  dnith's  program  for  the  gust  pertur¬ 
bation  case  and  are  referred  to  as  the  Whitehead-Smith  Theory. 

Using  the  theory  developed  by  Henderson  [  r> )  ,  Shen  [33]  obtained 
an  explicit  solution  for  the  unsteady  pressure  distribution  along  the 
chord  of  a  blade  in  a  two-dimensional  cascade  caused  by  spatial  inflow 
velocity  variations.  The  results  f ron  this  theory  are  referred  to  as 
the  Henderson-Shen  theory.  The  method  employed  is  a  Singularity  Vortex 
Method,  which  has  been  widely  used  in  thin-airfoil  theory.  However, 
the  solution  includes  the  effects  of  mean  incidence  angle,  blade  camber, 
and  space-chord  ratio,  which  are  important  parameters  in  the  analysis 
of  the  unsteady  lift  of  a  blade  in  a  cascade. 

Some  important  assumptions  of  their  model  are:  (1)  the  flow  is 
two-dimensional,  inviscid,  and  incompressible,  (2)  the  disturbance 

velocities  are  small  compared  with  the  mean  velocity,  W  ,  and  (3)  the 

m 

camber  of  the  foil  is  small  so  that  the  boundary  condition  is  con¬ 
sidered  to  be  on  the  chord  line  rather  than  the  camberline.  The 
resulting  unsteady  lift  L  can  be  expressed  as: 

L  =  C,  irpCA,  W  V  S 
L  f  m  x 


(7) 
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where  C  is  Che  complex  unsteady  life  coefficient  which  is  a  function 

L 

of  reduced  frequency  and  of  the  cascade  operation  parameters.  In  the 
present  experimental  study:  (1)  the  Mach  number  was  less  than  0.1 
for  all  test  cases;  (2)  the  designed  fundamental  harmonic,  ,  of  the 
distortion  screens  was  0.2;  and  (3)  the  maximum  camber-to-chord  ratio 
was  0.11  percent  at  the  midspan,  i.e.,  the  location  of  the  instrument 
blade  segment.  These  conditions  are  considered  appropriate  for  a 
linearized  analysis.  Also,  the  velocity  field  impinging  on  the  rotor 
was  designed  to  be  purely  axial,  and  the  one-inch  blade  segment  was 
centered  at  midspan  of  the  rotor  where  the  rlow  would  be  two-dimensional. 

3 . 3  Coupling  of  Unsteady  Lift  to  Acoustic  Modal  Pressures 

The  unsteady  lift  can  be  calculated  from  the  Henderson-Shen  predic¬ 
tion  program  as  a  function  of  spanwise  location.  For  plane  wave 
excitation,  the  velocity  distortion  harmonic  and  the  rotor  blade 
number  are  the  same,  so  the  lift  on  all  the  blades  will  be  in  phase. 

For  the  higher  order  (1,0)  mode,  there  is  a  phase  difference  between 
the  unsteady  lift  on  each  blade  since  the  distortion  harmonic  is 
different  from  the  blade  number.  However,  this  effect  is  considered 
in  the  lift  prediction  program  by  specifying  the  inter-blade  phase 
angle.  The  Tyler-Sofrin  analysis  is  then  employed  to  obtain  the 
acoustic  pressure  so  that  the  total  effect  can  be  obtained  sirapLy  by 
multiplying  the  lift  from  one  blade  by  the  total  number  of  blades. 

Due  to  the  blade  stagger  and  the  propagation  helix  angles,  the 
lift  vector  will,  in  general,  not  be  aligned  with  the  wave  front. 
Consequently,  the  lift  must  be  resolved  into  the  direction  of  the 
wave  front  as  shown  in  Figure  16.  Furthermore,  the  unsteady  lift  along 
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the  span  of  a  blade  must  be  weighted  spatially  by  the  radial  mode 

shape  function  E  (r)  and  summed  temporally  to  get  the  resultant 
tin 

unsteady  force.  This  resultant  unsteady  force  is  multiplied  by  the 
total  number  of  blades  and  divided  by  the  duct  area  to  obtain  the 
acoustic  modal  pressure  p  ,  i.e.  : 

rb 

r  E  (r)  L  sin  (£  +  a)  .  (8) 

*— 1  ran 
r 

a 

This  coupling  scheme  agrees  with  the  exact  solution  to  the  wave 
equation  as  described  by  Mani  [8);  the  comparison  is  described  in 
Appendix  E. 

The  flow  chart  of  the  noise  prediction  program  is  shown  in 

Figure  17.  In  the  experiment,  the  plane  wave  and  the  first  higher 

order  mode  (1,0)  were  of  main  interest.  For  the  nine-bladed  AFRF  rotor, 

the  velocity  distortion  harmonics  A„,  An,  and  A. _  excite  the  (+1,0) 

o  y  lu 

mode,  plane  wave,  and  (-1,0)  mode,  respectively.  Figures  18  and  19 
show  these  harmonics  versus  free-stream  velocity  at  three  radial 
locations.  These  values  were  read  into  the  prediction  program  and 
were  interpolated  to  obtain  the  vaLues  corresponding  to  the  operating 
conditions.  The  local  mean  axial  velocity  at  each  spanwise  segment 
was  also  obtained  from  interpolation  of  the  measured  velocity  profiles. 

Because  of  the  duct  termination  of  the  inlet  opening  and  the 
presence  of  rotor  blades,  a  standing  wave  will  result  from  the  net 
effect  of  the  outgoing  and  reflected  waves  between  these  two  planes. 

To  find  the  coupling  coefficient  between  the  unsteady  lift  and  acoustic 
modal  pressure,  it  is  necessary  to  decompose  the  measured  acoustic 
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Figure  17  Flow  Chart  of  Acoustic  Modal  Pressure  Preiction 


rier  Harmonics  of  Nine-Cycle  Screen  versus  Mean  Axial  Velocity 


Fourier  Harmonics  of  Ten-Cycle  Screen  versus  Mean  Axial  Velocity 
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pressures  at  BPF  into  their  modal  components  with  the  reflection  effects 
removed.  A  modal  decomposition  program  was  written  to  calculate  the 
modal  pressures  at  the  source  from  the  microphone  measurements.  A 
detailed  discussion  of  the  effect  of  reflection  and  modal  decomposition 
is  presented  in  Appendix  F. 

3 . 4  Comparison  of  Experimental  Results  with  Theoretical  Predictions 

3.4.1  Distortion  Screens 

The  desired  characteristics  of  the  flow  field  are  that  the  magni¬ 
tude  of  the  fundamental  harmonic,  non-dinensionalized  by  the  free-stream 
velocity,  has  a  value  of  0.2,  while  all  other  harmonics  are  relatively 
small.  In  order  to  obtain  a  strong  discrimination  between  the  plane 
wave  and  (-1,0)  mode,  the  fundamental  screen  circumferential  harmonic 
must  be  distinctly  larger  than  its  neighboring  ones.  Specifically,  Ag 
should  be  large  relative  to  Ag  and  A^g  for  the  nine-cycle  screen  and  A^g 
should  be  large  relative  to  Ag  and  Ag  for  the  ten-cycle  screen. 

The  results  of  the  Fourier  analysis  of  the  velocity  fields  down¬ 
stream  of  the  distortion  screens  are  shown  in  Figures  20  and  21.  It  was 
found  that  the  fundamental  harmonics  for  both  screens  were  near  0.2  as 
designed.  For  the  plane  wave  excitation,  nine-cycle  screen,  the  largest 
neighboring  distortion  harmonic  is  Ag ,  which  is  14.6  percent  of  Ag . 

Tne  (1,0)  mode  excited  by  Ag,  in  this  case,  will  decay  when  the  BPF  is 
below  the  cutoff  frequency  of  280  Hz.  Above  the  cutoff  frequency,  it 
will  have  a  contribution  to  the  acoustic  measurement. 

For  excitation  of  the  (-1,0)  mode,  the  most  significant  neighboring 
harmonic  is  Ag ,  which  is  36.3  percent  of  A^g  near  the  hub  and  13.8 
percent  near  the  tip.  The  plane  wave  excited  by  Ag ,  when  the  ten-cvcle 
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Figure  20  Amplitudes  of  Fourier  Harmonics  of  Nine-Cycle  Screen 


Amplitudes  of  Fourier  Harmonics  of  Ten-Cycle  Screen 
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screen  was  used  to  generate  the  (-1,0)  mode,  will  always  propagate 
and,  therefore,  its  presence  cannot  be  neglected.  The  presence  of  all 
these  modes  was  considered  in  the  decomposition  of  the  acoustic  pressures, 
i.e.,  the  measurements  were  decomposed  into  the  plane  wave  and  the 
(-1,0)  mode. 

The  largest  harmonic,  other  than  the  screen  fundamental,  appears 

near  the  hub.  For  the  nine-cycle  screen,  it  is  which  is  85  percent 

of  Aq  and  for  the  ten-cycle  screen,  it  is  A  which  is  55  percent  of 
y  o 

A^q.  However,  they  only  contribute  to  exciting  the  (+2,0)  mode.  In 
the  experiment,  BPF  was  near  280  Hz  which  is  well  below  the  cutoff 
frequency  of  the  (2,0)  mode,  i.e.,  550  Hz. 

Tyler  and  Sofrin  [10]  gave  an  expression  for  the  decay  rate  of 
the  acoustic  pressure  when  the  cutoff  ratio  is  less  than  one: 


A  plot  of  this  equation  for  the  various  modes  of  the  AFRF  is  shown 
in  Figure  22.  For  the  (2,0)  mode,  the  decay  rate  is  21  dB/ft  near 
280  Hz  and  each  higher  order  mode  decays  even  faster.  Therefore,  all 
the  modes  higher  than  the  (1,0)  mode  can  be  neglected  over  the  operating 
speed  range  investigated. 

Figures  20  and  21  also  show  that  the  worst  screen  performance 
occurred  near  the  hub.  However,  the  magnitude  of  the  unsteady  blade 
pressure  is  a  minimum  at  the  hub.  This  follows  because  the  blade  chord 
is  constant  along  the  span  so  the  reduced  frequency  is  about  the  same 
from  hub  to  tip  and,  thus,  the  section  lift  coefficient  is  about 


constant  along  the  span.  However,  the  unsteady  lift  is  proportional 


to  the  square  of  the  relative  velocity,  W  for  a  constant  ratio  of 

m 

W./W  and,  thus,  the  contribution  to  the  unsteady  pressure  near  the 
dm 

hub  is  small  because  the  local  relative  velocity  is  low  relative  to 
the  remainder  of  the  span. 

3.4.2  Unsteady  Lift 

Figures  23  to  25  show  typical  results  of  the  unsteady  lift 
measurement  for  the  signal  analyses:  spectrum  average  (function  4), 
waveform  average  (function  12),  and  waveform  average/FFT  (function  13) , 
respectively. 

In  Figure  25,  the  largest  harmonic,  N  =  16,  resulted  from  the 
mechanical  resonance  of  the  lift  gage  near  500  Hz,  since  1S54  rpm 
corresponds  to  494  Hz  for  the  sixteenth  harmonic  (1854  x  16/60  =»  494.4). 
This  peak  moved  to  higher  values  of  N  as  the  shaft  speed  decreased. 

Measurement  and  prediction  of  unsteady  lift  versus  rotor  speed 
for  different  test  conditions  are  presented  in  Figures  26  to  29.  It 
is  apparent  that  there  was  a  strong  back-reaction  influence  on  the 
unsteady  lift  magnitude  as  well  as  phase  for  both  the  nine-  and  ten- 
cycle  screens,  which  occurred  near  245  Hz  (1630  rpm  for  the  nine-cycle 
screen  and  1470  rpm  for  the  ten-cycle  screen).  The  theoretical  predic¬ 
tion  is  in  good  agreement  with  the  measurement  for  the  ten-cycle 
screen  except  for  the  variation  caused  by  the  back-reaction  effect. 
However,  the  discrepancy  between  measured  and  predicted  lift  values 
were  rather  large  for  the  nine-cycle  screen  cases. 

The  unsteady  lift  coefficients  at  different  spanwise  sections  for 
the  test  cases  are  listed  in  Table  III.  Operating  at  the  two  flow 
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TABLE  III 

UNSTEADY  LIFT  COEFFICIENTS  OF  THE  TEST  CASES 


Unsteady  lift  coefficient 


10.15 


Reduced 

Whitehead- 

Henderson-Shen 

Frequency 

Smith  theory 

<>  =  0.56 

A 

Mag.  Phase 

2.312 

0.099  130.96 

0. 1614  146.25 

2.601 

0.0854  144.95 

0.2780  128.05 

2.890 

0.1127  151.85 

0.1482  143.47 

2.473 

0.1744  126.82 

0.1633  128.44 

2.782 

0.1295  135.39 

0.3146  91.73 

3.091 

0.1587  147.33 

0.2020  133.40 

2.382 

0.2085  126.29 

0.1580  120.18 

2.679 

0.1695  129.66 

0.3663  85.85 

2.977 

0.1943  142.28 

0.2413  135.42 

2.217 

0.2362  127.03 

0.1695  111.49 

2.495 

0.2043  126.48 

0.4400  94.58 

2.772 

0.2225  137.67 

0.2592  141.51 

2.045 

0.2574  127.79 

0.2098  110.01 

2.300 

6.2311  12480 

0.5244  107.15 

2.556 

0.2447  134.09 

0.2769  148.81 

146.90 

130.03 

143.71 


130.20 

95.47 

133.76 


122.22 

89.08 

135.17 


142.29 

96.69 

140.46 


112.72 
103 . 24 
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coefficients  produced  different  mean  incidence  angles  and,  therefore, 
steady  lift  values.  The  steady  lift  coefficients  can  be  calculated 
by  using  the  correlation  of  lift  characteristic  for  blades  in  cascade, 
which  resulted  in  values  of  lift  coefficient  of  0.328  and  0.683  at 
midspan  for  <!>  =  0.70  and  0.56,  respectively.  The  mean  incidence 
angles  for  these  flow  coefficients  are  0.3  and  0.6  degrees,  respectively. 
Note  that  the  steady  loading  causes  a  decrease  in  the  magnitude  of 
the  unsteady  lift  coefficient,  C^,  by  about  ten  percent  using  the 
Henderson-Shen  theory.  Because  the  Whitehead-Smith  theory  does  not 
include  the  effect  of  steady  loading  on  the  blades,  the  value  of  C 

Lt 

remains  constant. 

When  the  nine-cycle  screen  was  installed  in  the  nine-bladed  fan, 
the  spacing  between  blades  was  equal  to  the  wavelength  of  the  distorted 
inflow,  i.e.,  the  interblade  frequency  approached  211,  which  is  referred 
to  as  an  "aerodynamic  resonance"  in  a  cascade  [5].  Henderson-Shen ' s 
prediction  scheme  predicts  a  significant  increase  in  C  near  resonance, 

Lt 

while  the  Whitehead-Smith  theory  does  not  have  a  change  in  at  this 
resonance  point.  The  measurements  showed  that  C  changed  in  magnitude 

Li 

and  phase  for  the  nine-cycle  screen  relative  to  that  of  the  ten-cycle 
screen,  and  both  analyses d  id  not  predict  this  reduction  or  change  in 
phase  angle. 

For  the  ten-cycle  cases,  the  agreement  between  the  Whitehead-Smith 
theory  is  excellent  when  the  rotor  speeds  were  less  than  1450  rpm. 

The  unsteady  lift  fluctuated  at  higher  rotor  speeds,  which  apparently 
is  the  result  of  the  back-reaction  of  the  strong  acoustic  field. 

The  Henderson-Shen  theory  overestimates  by  15  <115  in  the  magnitude; 


however,  the  prediction  of  the  phase  angle  is  good.  The  larger  discrep- 
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ancies  in  the  phase  angle  from  the  Whitehead-Smith  theory  for  the 
cases  <i>  =  0.56  are  expected,  since  the  mean  incidence  angle  was  6.0 
degrees,  which  departed  from  the  zero  steady  loading  of  their  model. 

Bruce  [12]  has  conducted  a  series  of  experiments  in  the  AFRF  to 
measure  the  unsteady  lift  and  moment  generated  by  the  disturbance  flow 
field.  The  comparison  of  the  data  measured  by  Bruce  and  the  predictions 
[12,  27]  also  showed  the  theories  tend  to  overestimate  the  magnitude  of 
the  unsteady  lift,  while  the  measured  and  predicted  phase  angles  were 
in  good  agreement  even  near  the  aerodynamic  resonance.  However,  the 
reduced  frequencies  in  Bruce's  experiment  were  0.2  to  2.0,  which  were 
lower  than  the  range  of  reduced  frequencies,  2.7  to  3.0  used  in  this 
study.  In  particular,  the  reduced  frequency  near  the  aerodynamic 
resonance  was  1.3  in  Bruce's  experiment  and  2.7  in  this  experiment . 

3.4.3  Acoustic  Pressure 

Figure  30  shows  the  sound  pressure  level  measured  at  microphone 
position  one  over  the  speed  range  for  the  case  <J>  =  0.70.  The  measure¬ 
ments  show  that  the  acoustic  pressure  measured  with  the  nine-cycle 
screen  is  practically  the  same  over  a  large  range  of  rotor  speeds. 

For  the  ten-cycle  screen,  a  cutoff  phenomenon  was  observed  near  1850  rpm 
as  predicted.  Thus,  the  acoustic  pressure  was  dominated  by  the  plane 
wave  mode  for  the  nine-cycle  screen  and  by  the  (-1,0)  mode  for  the 
ten-cycle  screen,  as  designed. 

The  complete  data  of  the  magnitude  and  phase  angle  of  the  acoustic 
pressures  measured  at  four  microphone  locations  are  listed  in  Table  IV 
and  V.  These  data  were  used  as  the  input  to  the  modal  decomposition 


program  described  in  Appendix  F.  The  resulting  decomposed  modal 


TABLE  IV 


SOUND  PRESSURE  LEVEL  MEASUREMENT  —  TEN-CYCLE  DISTORTION  SCREEN 
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pressures  are  shown  in  Figures  31  to  34.  The  plane  wave  was  found 
to  be  larger  than  the  (-1,0)  mode  by  6  to  11  dB  for  the  nine-cycle 
screen,  and  the  (-1,0)  mode  was  larger  than  the  plane  wave  by  16  to  17 
dB  for  the  cases  above  the  cutoff  frequency.  The  discrimination  between 
the  two  modes  had  a  low  value  at  cutoff,  about  6  dB  for  the  nine-cycle 
screen  and  4  dB  for  the  ten-cycle  screen. 

The  reason  for  the  poor  discrimination  between  the  two  modes  at 
cutoff  may  be  caused  by  the  reflection  coefficient  of  the  rotor  being 
large  which  may  be  associated  with  the  observed  back-reaction  effect. 

The  reflection  coefficient  of  the  rotor  at  all  other  conditions  is 
theoretically  very  low  and,  for  that  reason,  it  was  not  included  in  the 
decomposition  scheme. 

It  is  important  to  know  whether  the  SPL  measured  by  the  microphone 
was  influenced  by  the  unsteady  pressure  of  the  aerodynamic  field  or  not. 
This  can  be  checked  by  considering  the  unsteady  aerodynamic  pressure  on 
the  blades  as  a  potential  field  of  a  flow  passing  a  wave-shaped  wall. 
Shapiro  [34]  indicated  that  the  decay  rate  of  hydrodynamic  pressure 
of  this  field  in  the  direction  perpendicular  to  the  wall  is 
exp[-2~>/l-M^  x/S],  where  M  is  the  Mach  number  of  the  flow,  and  S  is 
the  wavelength  of  the  wall,  which  corresponds  to  the  spacing  between 
the  blades  in  this  application.  Setting  S  =  0.625  ft  at  the  tip, 

M  =  0.17  at  2100  rpm,  the  highest  circumferential  speed  of  the  test 
cases,  the  resultant  decay  rate  is  -86.0  dB/ft.  The  highest  value  of 
unsteady  aerodynamic  pressure  on  the  outer  one-inch  of  the  blade  span 
is  150  dB  at  2100  rpm.  Therefore,  the  above  formula  indicates  it  will 
decay  to  -7.7  dB  when  it  reaches  the  nearest  microphone  position  of  22 
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Figure  32  Decomposed  Modal  Pressure  of  Nine-Cycle  Screen 
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Figure  34  Decomposed  Modal  Pressure  of  Ten-Cycle  Screen  (}  -  0.56) 


73 

inches  from  the  rotor.  This  value  is  well  below  the  typical  measured 
SPL  value  of  130  dB.  Thus,  there  will  be  no  contribution  of  the 
unsteady  hydrodynamic  field  to  the  acoustic  measurement. 

The  values  of  the  reflection  coefficient  of  the  plane  wave  at 
the  duct  opening,  R,  obtained  from  the  decomposition  program  are  plotted 
in  Figure  35.  It  can  be  seen  that  R  is  always  less  than  unity  as 
expected;  however,  it  varies  with  the  distortion  screens  and  flow 
coefficient.  The  solid  line  represents  the  prediction  of  R  by  the 
Morfey  theory  as  described  in  Appendix  F.  The  large  discrepancies 
between  measurement  and  prediction  indicates  that  the  reflection 
mechanism  could  be  more  complicated  than  the  theoretical  model  considered 
in  the  present  analysis. 

The  prediction  of  the  modal  sound  pressure  level  (SPL)  for  different 
values  of  unsteady  lift  coefficient,  C^,  is  also  plotted  in  Figures  31 
to  34.  For  the  plane  wave  excitation,  the  predicted  SPL  was  found  to 
be  always  higher  than  the  decomposed  value.  The  discrepancy  was: 
about  15  dB  using  from  the  Henderson-Shen  theory,  and  about  9  dB 
using  from  the  Whitehead-Smith  theory.  Prediction  of  the  plane  wave 
pressure  using  the  unsteady  lift  coefficients  obtained  from  the  previous 
lift  gage  measurement  was  also  computed,  and  it  was  still  higher  than 
the  decomposed  SPL  by  7  dB.  For  the  (-1,0)  mode  excitation,  the 
decomposed  modal  pressure  can  be  obtained  only  when  the  rotor  speed  was 
above  cutoff,  since  the  (-1,0)  mode  decays  quickly  along  the  duct  when 
the  cutoff  ratio  is  less  than  unity.  The  predicted  SPL  using  the 
derived  in  the  three  above  ways  underestimated  the  modal  pressure  for 
most  cases.  The  agreement  of  the  phase  angle  between  prediction  and 


measurement  was  not  good  for  both  screens. 
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It  is  difficult  to  explain  the  reasons  for  the  discrepancy  between 
the  predictions  of  the  acoustic  pressure  and  the  unsteady  lift.  Possibly, 
the  reasons  are  that:  (1)  The  measured  SPL  is  a  standing  wave  which 
results  from  the  acoustic  reflections,  and  varies  from  (1-R)  to  (1+R) 
times  the  source  pressure  depending  on  the  measurement  location.  The 
large  fluctuations  of  reflection  coefficient  in  Figure  35  shows  that 
a  more  detailed  measurement  of  the  acoustic  field  in  the  duct  is 
necessary  to  get  a  reliable  decomposed  modal  pressure  at  the  source. 

(2)  The  effect  of  the  boundary  layer  and  the  tip  clearance  might  not 
be  negligible  in  calculating  the  overall  SPL  at  BPF.  The  local  velocity 
is  much  smaller  than  the  mean  axial  velocity  in  the  boundary  layer  ne,  r 
the  tip  and  hub  of  the  rotor,  consequently,  the  unsteady  lift  would  be 
lower  because  the  lift  is  proportional  to  the  square  of  local  axial 
velocity.  The  contribution  to  unsteady  lift  near  the  tip  is  important, 
since  the  linear  circumferential  velocity  is  a  maximum  at  that  point. 
Therefore,  the  effect  of  tip  clearance  and  the  velocity  discrepancy 
might  have  a  strong  influence  when  predicting  the  overall  SPI.  in  the 
duct.  For  example,  the  percentage  of  the  total  acoustic  pressure  that 
is  associated  with  the  outer  one  inch  of  the  blade  span  is  67%.  This 
is  roughly  a  factor  of  -3.5  dB.  This  would  reduce  the  predicted  SPL 
by  this  same  amount  and  improve  the  agreement  between  all  the  predic¬ 
tions  and  measurements  for  the  plane  wave  case.  However,  for  the 
spinning  mode  case,  the  measured  acoustic  values  are  distinctly  greater 
than  the  predicted  values.  Furthermore,  the  measured  and  predicted 
phase  angle  is  also  in  distinct  disagreement.  Again,  because  the 
decomposed  mode  pres:  ire  depends  upon  the  reflection  coefficients  and 
the  theoretically  assumed  distribution  of  the  node  pressure  in  the 
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duct,  the  possibility  of  error  exists.  Acoustic  measurements  at 
more  locations  along  the  duct  would  resolve  the  question  of  whether 
the  acoustic  behavior  is  as  predicted  in  Appendix  F  or  whether  there 
is  an  error  in  the  assumptions  of  the  analysis  used  to  predict  the 
acoustic  modal  pressure  from  the  unsteady  lift  values.  Because  there 
is  observed  a  strong  back-reaction  of  the  acoustic  pressure  on  the 
unsteady  lift  values,  it  is  quite  possible  that  some  modification  to 
the  theory  will  be  necessary. 


CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

4 . 1  Conclusions 

The  following  conclusions  are  drawn  from  this  experimental  investi¬ 
gation  : 

1.  In  the  comprehensive  experiment  described,  velocity  perturba¬ 
tions,  fluctuating  blade  lift,  and  overall  noise  levels  were 
measured  in  order  to  produce  a  better  understanding  of  the 
basic  fan  noise  mechanism.  By  using  distortion  screens,  a 
plane  wave  and  a  spinning  wave  were  excited  as  designed. 

There  was  sufficient  signal-to-uoise  ratio  for  the  unsteady 
lift  and  acoustic  pressure  measurements ;  thus,  it  was  con¬ 
firmed  that  the  AFRF  is  a  suitable  facility  for  studying 
acoustic  duct  modes. 

2.  Distortion  screens  to  produce  +20  percent  sinusoidal 
velocity  profiles  were  designed,  built,  and  tested.  Good 
agreement  was  achieved  between  design  and  measurement. 

3.  The  lift  gage  in  the  instrumented  rotor  blade  was  calibrated 
as  described  in  Appendix  G,  and  its  resonance  frequency  was 
found  to  be  above  the  frequency  range  of  interest  in  this 
study.  The  electrical  noise  from  the  slip  ring  was  insigni¬ 


ficant  compared  with  the  unsteady  lift  signals. 
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4.  Using  a  phase-locked  ensemble  averaging  technique,  the 
amplitude  and  phase  angle  of  the  unsteady  lift  and  acoustic 
pressure  were  measured  with  good  repeatability. 

5.  A  cutoff  phenomenon  was  observed  for  the  (-1,0)  mode.  Most 
of  all,  a  significant  back -reaction  was  observed  near 
cutoff  on  the  aerodynamic  lift  for  both  the  plane  wave  and 
(-1,0)  mode  excitation,  which  had  been  suggested  by  Doak 
[23]  and  for  which  no  experimental  confirmation  has  been 
reported. 

6.  A  noise  prediction  program  and  a  modal  decomposition  program 
were  written  which  gave  a  reasonable  trend  for  the  amplitude 
of  the  decomposed  modes.  However,  the  agreement  with  the 
theoretical  prediction  was  poor.  The  discrepancies  are 
believed  to  have  resulted  from  the  reflection  model  used  in 
this  study  not  being  sufficiently  refined  to  describe  the 
resulting  acoustic  standing  wave. 

7.  The  Whitehead-Smith  and  Henderson-Shen  theories  were  success¬ 
ful  in  predicting  the  unsteady  lift  on  the  blade  at  midspan 
for  the  ten-cycle  screen.  However,  both  models  did  not  provide 
a  satisfactory  prediction  for  the  case  of  "aerodynamic 
resonance"  when  the  nine-cycle  screen  was  installed. 

8.  The  modal  coupling  scheme  used  in  the  prediction  of  SPL  is 

direct,  i.e.,  there  is  no  amplification  in  amplitude  and 
phase  shift  from  unsteady  aerodynamic  lift  to  acoustic  modal 
pressure.  Since  the  result  of  modal  pressure  from  the  decompo¬ 
sition  program  is  doubtful,  the  coupling  coefficient  in  both 
magnitude  and  phase  is  still  undeterm ined .  ^ 


79 


4 . 2  Recommendat ions 

The  recommendations  for  future  work  are: 

1.  A  fuller  knowledge  of  the  effect  of  acoustic  reflection  is 
necessary  to  decompose  the  measured  acoustic  pressures  and 

to  couple  the  unsteady  lift  to  acoustic  duct  modes.  By  using 
a  larger  number  of  duct  positions  (both  axially  and  circum¬ 
ferentially)  to  measure  the  amplitude  and  phase  of  sound 
pressure,  the  modal  structure  [21,  34] ,  the  reflection 
coefficient,  and  the  distance  between  successive  peaks  [6,  11] 
can  easily  be  obtained  with  greater  accuracy. 

2.  A  detailed  measurement  of  unsteady  pressures  on  the  rotor 
blade  by  using  miniature  pressure  transducers  on  both  the 
pressure  and  suction  sides  should  provide  more  information 
on  the  back-reaction  of  the  acoustic  pressure  on  the  aero¬ 
dynamic  field. 

3.  This  investigation  used  a  ten-cycle  screen  which  produced  a 
(-1,0)  rotating  mode.  An  eight-cycle  screen  should  be 
similarly  tested  as  it  will  produce  a  (+1,0)  mode.  The 
difference  in  direction  of  mode  rotation  will  change  the 
reflection  coefficient  of  the  rotor  and  may  provide  some 
additional  insight  of  the  back-reaction  behavior  near  cutoff. 

4.  An  experimental  study  of  the  cutoff  behavior  and  back-reaction 
of  higher  order  modes  other  than  (1,0)  mode  is  desirable.  By 
increasing  the  number  of  rotor  blades  and  stagger  angle,  it 
would  be  possible  to  reach  the  cutoff  frequencies  of  these 
higher  order  modes  and  still  maintain  a  nominal  zero  incidence 
condition  within  the  capacity  of  the  AKRF. 
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To  confirm  the  validity  of  the  theoretical  prediction  of 
unsteady  lift  coefficient  other  than  the  "aerodynamic 
resonance"  condition,  a  wider  range  of  reduced  frequencies 
should  be  tested.  This  can  be  obtained  by  changing  the 
stagger  angle  and  the  distortion  screens.  For  a  detailed 
check  of  the  validity  of  the  lift  theories  near  aerodynamic 
resonance,  it  is  also  necessary  to  select  a  wider  range  of 
blade  numbers  with  the  same  number  of  cycles  in  the  distor¬ 
tion  screen  to  measure  the  unsteady  lift  response. 

To  confirm  that  the  hydrodynamic  field  does  indeed  decay  to 
a  sufficiently  low  level  at  the  microphone  measurement 
locations,  more  measurements  close  to  the  rotor  are  needed. 

A  theoretical  model  which  can  describe  the  character ' sties 
of  the  back-reaction  of  the  acoustic  field  on  the  aerodynamic 
field  should  be  developed. 
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APPENDIX  A 


DESIGN  OF  THE  ROTOR 

The  steady  state  design  of  the  AFRF  6-inch  chord  rotor  was 
conducted  by  using  Lieblien's  [36]  correlations,  which  are  based  on 
a  minimum  loss  condition.  At  this  condition,  the  actual  total  pressure 
ratio  across  the  stage  is  maximized.  The  rotor  blade  shapes  were 
specified  as  having  a  circular-arc  camberline  with  a  10  percent  maximum 
thickness-to-chord  ratio. 

The  approach  used  was  to  apply  a  correction  for  spanwise  effects 
to  a  mean  radius  design.  The  mean  radius  design  consists  of  four 
aspects:  (1)  velocity  diagram,  (2)  number  of  rotor  blades,  (3)  the 

blade  geometry  (camber  and  stagger  angle),  and  (4)  the  stall  limit. 

A. 1  Mean  Radius  Design 

A.  1.1  Velocity  Diagram 

The  velocity  triangle  at  the  rotor  inlet  is  constructed  using  the 
axial  velocity  calcualted  from  the  continuity  equation  for  the  design 
mass  flow  rate,  and  the  blade  rotational  velocity  U,  determined  from  the 
design  shaft  speed.  The  change  of  the  absolute  circumferential  velocity 
can  be  calculated  from  Euler's  pump  equation  for  a  specified  energy 
input  from  the  drive  motor  at  the  design  condition;  this  permits 
construction  of  the  velocity  triangle  at  the  rotor  outlet. 
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A. 1.2  Number  of  Rotor  Blades 

Howell  [37]  correlated  the  space-to-chord  ratio,  S/C,  and  fluid  angle 
at  "nominal  deflection"  (80  percent  of  the  deflection  angle  at  which 
the  blades  start  to  stall)  as: 


1-55  ,  pi  x  1 

tan  3  -  tan  8  '  1.5  ’ 

'  i  e  y 


where  6  and  3  are  the  fluid  angles  relative  to  the  rotor  at  inlet  and 
i  e 

outlet,  respectively,  that  have  been  determined  from  the  velocity 
triangle.  Since  B  =  2:^/5  is  not  an  integer,  B  was  taken  as  the 
closest  integer. 


A.  1 . 3  Blade  Geometry 

The  camber  angle  correlated  by  Lieblien  [36]  can  be  expressed  as: 

o  -  (3i-V  -Ksh^  “Vio-  W 

1  -  m  +  n  ’  (‘A-; 

where  m,n  are  correlation  constants,  and  and  K.  are  shape  correction 

factors.  The  incident  angle  i  and  deviation  angle  6  can  be  obtained 


i  =  i  +  no  and  6  =»  6  +  ma  .  ( 

o  o’ 

where  iQ,  6  ,  m,  and  n  are  plotted  as  correlation  curves.  By  defini¬ 
tion,  the  blade  angles  8!  and  3'  are: 

l  e 


8!  =*  3.  -  i  and  B'=8-5 
li  e  e 
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and  the  stagger  angle  is: 


3!  +  £' 


(A4) 


A . 1 . 4  Stall  Limit 

Lieblien  defined  the  diffusion  factor,  D,  as 


D  =■  1 


* 

cos 

6i 

's' 

cos  6 

1 

cos 

0 

e 

+ 

Cj 

2 

[tan  8.  -  tan  8  ] 
i  e 


For  an  unstalled  condition,  it  is  required  that  D  <  0.6. 


(A5) 


A. 2  Spanwise  Variation 

The  hub-tip  radius  ratio  for  the  AFRF  is  0.442,  which  is  lower  than 
the  limit,  0.7,  for  a  two-dimensional  approach,  so  the  radial  variation 
of  flow  also  should  be  considered.  To  avoid  swirl  in  the  flow,  a  free- 
vortex  design  was  employed,  i.e.: 


and 


Vr  r  =  vr  rm  “  0  at  rotor  inlet 

i  i 

Vr  r  =  Vr  rm  =  const.  at  rotor  exit 

e  e 


The  fluid  angles  at  different  radii  can  be  calculated  from 


and 


3^  =*  tan  1  (r-o/V^) 

B  =■  tan  1  [(rui  -  V  )/V  1 
e  r  x 

e 

a  =■  tan  1  (V  /V  ) 
a  r  x 


(A6) 


Finally,  the  blade  angle  at  each  radius  can  be  determined  in  the  same 
way  as  the  midspan  analysis,  and  unstalled  conditions  checked  at  each 

radius . 
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APPENDIX  B 


SELECTION  OF  MICROPHONE  LOCATIONS 

The  objective  of  the  microphone  measurement  is  to  decompose  the 
measured  acoustic  pressure  into  duct  modal  components  and  to  confirm 
that  the  generated  pattern  is  the  (0,0)  or  (-1,0)  mode  as  designed. 

The  reflection  effect  of  rotor  and  duct  inlet  opening  was  checked 
by  the  reflection  program  described  in  Appendix  F,  and  it  was  found 
that  a  steady  state  pressure  distribution  is  calculated  after  the 
first  reflection  by  the  duct  inlet  opening  for  the  (0,0)  and  (-1,0) 
modes  for  all  conditions  including  the  more  critical  condition  near 
cutoff  in  the  AFRF  as  shown  in  Figure  F2.  Thus,  only  the  first 
reflection  is  considered  in  the  reverse  problem  of  acoustic  decompo¬ 
sition. 

To  determine  the  number  of  microphone  measurement  positions 
needed  for  modal  decomposition,  consider  two  cases:  (i)  the  cutoff 
ratio  is  greater  than  unity,  and  the  (-1,0)  mode  will  propagate,  so 
there  are  four  unknowns:  p^,  p^  [complex  pressure  for  (0,0)  and 
(-1,0)  modes],  and  Rq  and  R^  [complex  refraction  coefficients  of  the 
inlet  opening  for  either  mode].  Therefore,  at  least  four  measurement 
points  are  needed  to  decompose  the  modal  pressure  p^,  p^  and  compare 
Rq,  R^  with  the  theoretical  prediction.  In  case  (ii),  the  cutoff  ratio 
Y  <  1,  only  the  plane  wave  propagates,  so  p^  and  are  eliminated  and 
only  two  measurement  points  are  needed. 


To  discriminate  between  the  various  modes,  the  microphone  axial 
location  should  be  optimized.  This  is  accomplished  by  locating  the 
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microphone  where  the  desired  mode  has  a  higher  pressure  magnitude 
relative  to  all  other  modal  components.  In  the  test  program,  only 
(0,0)  and  (-1,0)  modes  are  of  interest.  The  plots  of  | pQ |  -  | p^ | 
versus  axial  position  were  run  for  some  test  cases.  The  optimized 
measurement  location  for  the  plane  wave  is  where  the  quantity  ( I PQ I 
|p^|)  is  a  maximum,  and  for  the  (-1,0)  mode,  where  the  quantity 
|pQ  |  -  |p1  I)  is  a  minimum.  However,  the  standing  wave  pattern  is  a 
function  of  frequency,  i.e.,  the  optimized  locations  vary  with  y ,  so 
the  four  locations  were  selected  for  each  mode  near  the  optimized 
positions  over  the  cutoff  frequency  ratios  used  in  the  test.  Microphone 
locations  close  to  the  rotor,  which  are  in  the  near  field  (<_1.0  ft) 
were  avoided  to  ensure  that  the  near-field  hydrodynamic  pressure  had 
decayed  sufficiently  so  that  it  did  not  contribute  to  the  acoustic 
measurement . 

At  the  same  axial  location,  different  circumferential  locations 
introduce  a  phase  shift  in  the  pressure  measurement  for  the  higher 
order  spinning  mode  when  y  >  1.  There  is  no  circumferential  phase 
dependence  for  the  plane  wave.  Consequently,  the  four  selected  axial 
locations  were  staggered  circumferentially  to  distinguish  one  mode  from 
the  other  by  examining  the  phase  difference  between  each  measurement 


location . 


APPENDIX  C 


DESIGN  OF  DISTORTION  SCREENS 

Tyler  and  Sofrin  [10]  indicated  that  the  circumferential  order  of 
an  acoustic  mode  generated  by  rotor-stator  interaction  can  be  expressed 
as : 

m  =  nB  +  kV  ,  (Cl) 

where 

B  =  number  of  rotor  blades, 

V  =  number  of  stator  blades  or  distortions  at  rotor  inlet  flow, 
n  =  BPF  harmonics 

and 

k  =  ...  -1,  0,  1  ...  arbitrary  integer. 

For  the  nine-bladed  rotor  in  the  AFRF,  this  equation  shows  that 
the  plane  wave  (0,0)modeand  the  first  higher  order  (-1,0)  mode  can 
be  generated  by  a  nine-  and  a  ten-cycle  distortion,  respectively. 

Nine-  and  ten-cycle  distortion  screens  were  employed  in  the 
experiment.  The  screens  were  designed  to  produce  sinusoidal  velocity 
profiles  with  ±20  percent  variation  in  the  mean  axial  velocity  V^_. 
Sinusoidal  velocity  distortions  are  desired  because  the  aero-acoustic 
analysis  is  linear.  Consequently,  knowing  the  behavior  for  a  pure  sine 
wave  permits  prediction  for  more  complicated  waveforms  using  superposition. 
A  distortion  magnitude  of  20  percent  was  selected  as  it  is  a  significant 
distortion,  but  still  small  enough  to  apply  linearized  theories  in 
predicting  the  response  of  the  rotor. 
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For  the  rotor  inlet  flow,  the  spatially-fixed  axial  velocity 
can  be  decomposed  into  a  Fourier  series  as: 


x  AG 


E  ^ 


i(nG  +  <j>  ) 


The  desired  velocity  profile  should  be  dominated  by  either  the 

ninth  or  the  tenth  harmonic  in  the  experiment. 

On  passing  through  a  screen,  the  total  pressure  of  the  flow  is 

reduced.  This  pressure  drop  A?  can  be  described  by  a  dimensionless 

resistance  coefficient  K,  which  is  based  on  the  local  velocity  V 

J  o 

normal  to  the  screen: 


1  v2 

y  P  v. 
i  o 


McCarthy  [38]  developed  an  analytical  solution  for  a  steady 
moderately  sheared  three-dimensional  flow  past  wire  screens.  The 
solution  expressed  the  velocity  downstream  of  the  screen  as  a 
function  of  flow  resistance,  K.  Substituting  the  equation  into  the 
Fourier  series  leads  to  the  expression 


1.02(1  +  X)  1  +  6X 

A  stn  nO  =  - TTT\  ^ - y— 

f  (1  +  X3)273  °  6X^ 


where 


1/2 

X  =  (1  +  K)  ' 


X  =  (1  +  K  .  ) 
o  m  in 


1  +  6X  3  A.( 1  +  X  3) 

_ t  o 

Yo  -  2  1.02(1  +  X  ) 

OX  o 

o 


2 
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The  resistance  coefficients  of  the  support  screens  K  .  are 

mm 

0.343  and  0.236  for  the  nine- and  the  ten-cycle  screens,  respectively. 
Setting  the  magnitude  of  the  fundamental  harmonic  =  0.2,  the  solution 
of  this  equation  gives  the  desired  variation  of  K  with  0  as  shown  in 
Figures  Cl  and  C2. 

Various  segments  of  screens  with  different  resistance  coefficients 
were  overlaid  to  obtain  an  approximation  of  the  calculated  resistance 
variation.  The  approximation  was  determined  by  equally  dividing  the 
difference  in  area  between  the  step  changes  in  resistance  and  the 
desired  variation.  Ideally,  the  more  step  changes  used  in  the  K  distri¬ 
bution,  the  better  will  be  the  approximation  to  a  sinusoidal  velocity 
distribution.  However,  due  to  the  limitation  of  the  available  screen 
stock  and  the  difficulty  in  cutting  large  spacing  screens  into  small 
segments,  only  three-  and  one-step  changes  of  K  were  chosen  in  the  design 
of  the  screens.  A  more  detailed  design  procedure  was  presented  by 


Bruce  [39]. 
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Figure  Cl  Theoretical  and  Designed  Variation  of 

Resistance  Coefficient  for  Nine-Cycle  Screen 
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Figure  C2 


Theoretical  and  Designed  Variation  o£ 
Resistance  Coefficient  for  Ten-Cycle  Screen 


APPENDIX  D 

SIGNAL  PROCESSING  CONSIDERATIONS 

The  frequency  range  selected  on  the  analyzer  was  500  Hz,  since 
the  value  of  the  BPF  over  the  test  range  was  less  than  310  Hz  for  all 
test  conditions.  The  sampling  period  at  this  frequency  range  setting 
was  100  msec,  which  was  a  "time  duration  window"  long  enough  to  cover 
the  lowest  rotor  speed  condition,  i.e.,  742  rpm  or  1.35  msec/rev. 

The  aliasing  filter  was  disconnected  in  the  phase-locked  average 
tests,  since  a  phase  shift  would  be  introduced  by  the  filter  as  a 
function  of  input  frequency  that  was  undesirable.  Because  the  aliasing 
frequencies  were  higher  than  400  Hz,  which  is  beyond  the  frequency 
range  of  interest  for  the  test,  no  error  was  introduced  from  aliasing 
of  the  spectral  analyses. 

The  frame  size  of  the  signal  processor  is  1024  for  the  calculation 
of  the  spectrum  by  the  FFT  in  function  13.  These  data  points  were 
triggered  by  the  tracking  ratio  tuner  with  a  sampling  frequency  of 
1024/N  times  the  rotor  speed,  where  N  is  an  integer.  When  N  -  1,  the 
magnitude  and  the  phase  are  transformed  by  the  averaged  waveform  exper ienced 
by  the  instrumented  rotor  blade  in  one  revolution.  As  N  increases,  the 
sampling  rate  decreases  and  the  averaged  waveform  shrinks  in  the  time 
domain;  thus,  the  transformed  harmonics  of  the  shaft  speed  will  expand 
along  the  frequency  axis,  giving  a  better  resolution  to  distinguish  the 
fundamental  harmonic  and  the  neighboring  ones.  However,  although  there 
are  N  more  revolutions  in  the  1024  display  than  when  N  =  1,  the  accuracy 


-  - 
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is  decreased  since  there  are  fewer  sampling  points  per  revolution. 
Figure  1)1  shows  the  effect  of  changes  in  N  in  the  transformed  plane 
of  the  unsteady  lift  measurement.  The  value  N  =  4  was  taken  for  the 
lift  gage  measurements  to  obtain  better  resolution  when  recording  the 
data. 

Note  that,  although  better  resolution  of  the  harmonic  of  interest 
is  obtained,  the  level  of  the  harmonic  is  not  affected  by  the  N  setting. 
Since  the  acoustic  pressures  were  dominated  by  a  pure  tone  at  BPF, 
resolution  was  not  important  to  separate  tbe  harmonics;  thus,  N  =  1 
was  taken  in  the  acoustic  measurement. 


APPENDIX  E 


COMPARISON  OF  NOISE  PREDICTION  SCHEME  AND  EXACT  SOLUTION  TO 
SOUND  GENERATION  EQUATION 


Lighchill  [40]  firsC  derived  a  general  linearized  equation  for 
sound  generation  as: 


2 

c~ 


ox .  3:< . 

i  3 


(El) 


The  left-hand  side  of  the  equation  is  the  wave  equation,  and  the 
right-hand  side  can  be  regarded  as  a  collection  of  acoustic  source 
terms,  where  0,  the  first  term  on  the  right-hand  side,  is  the  rate  of 
introduction  of  mass  per  unit  volume  and  represents  the  sound  generated 
by  the  volume  displacement  effects  of  the  blades,  analogous  to  an 
acoustic  tr.oncpole.  F  is  the  external  force  per  unit  volume  acting  on 
the  fluid;  thus,  the  second  term  can  be  interpreted  as  a  dipole 
source  from  the  fluctuating  forces  exerted  on  the  rotor.  T,j  is  the 
stress  tensor  consisting  of  Reynolds  stress  and  fluid  pressures. 

The  last  term  represents  a  quadrupole  source  arising  from  turbulence 
in  the  flow  field. 

In  the  generation  of  fan  noise  at  low  Mach  number,  the  signficant 

2 

source  term  Ls  the  fluctuating  force  F^.  By  putting  p  =  c  o',  the 
equation  can  be  easily  converted  to  pressure,  which  is  applied  in  its 
most  practical  form  as: 


(E2) 


Mani  [8]  considered  a  temporal  harmonic  pressure  wave  and 
replaced  the  unsteady  force  on  the  cascade  by  a  row  of  delta  functions 
in  cylindrical  coordinates.  The  solution  of  Equation  (E2)  was  found 
to  be: 


p  =  )  p  E  (r)  exp  it ~  nit  i  U10  +  k  x  +  A  ] 

L  ran  ran  mn  ran 

ran 


where 


p  =  2tt 
mn 


VF  *  E  (r)  rdr 
mn 


and  VF  =  [ cos (£)  •  <S’(x)  +  i  sin(C)  ^  6(x)]  .  (E 3) 

In  the  expressions  of  Lighthill  and  Mani,  F  is  a  force  per  unit 
volume.  In  the  noise  prediction  program  used  in  this  thesis,  the 
force  is  the  lift  force  per  unit  span.  In  Mani’s  formulation,  the 
spatial  divergence  of  the  force  per  unit  volume  was  obtained,  resolved 
in  the  direction  of  the  wavefront,  weighted  by  the  radial  distribution 
of  mode  pressure,  and  integrated  over  the  duct.  In  the  noise  prediction 
scheme  used  in  this  thesis,  the  unsteady  lift  per  unit  span  was  resolved 
in  the  direction  of  the  wavefront,  integrated  along  the  span  and  weighted 
by  the  radial  mode  shape,  and  divided  by  the  duct  area.  Composing  the 
dimensions  of  the  quantities  used  in  each  prediction  shows  that  the 
result  is  a  unit  of  pressure. 


Mani : 
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Thesis : 


*  ' 

F 

x  E  x 

■  • 

L 

.  L  . 

mn 

. 

Note  the  unsteady  lift  force  is  assumed  to  couple  directly  in  magnitude 
with  no  temporal  phase  shift  to  the  acoustic  wavefront. 


APPENDIX  F 


THE  EFFECT  OF  REFLECTION  AND  ACOUSTIC  MODE  DECOMPOSITION 

The  discontinuity  of  the  duct  at  the  inlet  opening  and  the  rotor 
blade  causes  sound  waves  to  be  reflected  between  these  two  faces.  The 
measured  sound  pressure  at  a  fixed  position  is  then  the  net  effect  of 
the  incoming  wave  generated  at  the  source  and  the  reflected  waves 
which  travel  back  and  forth  between  the  duct  opening  and  the  rotor 
blade  row. 

The  pressure  pattern  of  higher  order  modes  (m,0)  can  be  written 

as : 


(i) 

below  cutoff 

P  = 

P  (r) 
ao 

-kx  i  (cot  -  0) 
e  e 

> 

(ii) 

at  cutoff 

P  = 

P  (r) 
mo 

ei(wt  -  0) 

J 

(iii) 

above  cutoff 

p  = 

P  (r) 
mo 

ei(wt  -  akx  -  0) 

• 

(FI) 

Examining  these  expressions,  when  the  frequency  is  below  or  at 
cutoff,  the  wave  front  is  in  the  0  direction,  which  is  parallel  to  the 
duct  opening,  so  there  will  be  reflection  in  these  cases.  For  the 
plane  wave  which  always  travels  axially  and  spinning  higher  order  modes 
above  their  cutoff  frequencies  which  travel  in  a  helix  path,  the  effect 
of  reflection  must  be  considered- 


F. I  Reflection  of  the  Duct  Opening 
Consider  an  oblique  incident  wave  with  helix  angle 
written  as: 

a  p  jUt  -  kd  ) 


it  can  be 


( F2 ) 


103 


where  cl  is  the  wave  front  vector,  d  =  x  cos  a  +  y  sin  a,  and  P  is  the 
complex  amplitude  of  the  incident  wave. 

The  angle  of  incidence  must  be  equal  to  the  angle  of  reflection, 
so  the  reflected  wave  is  written  as: 

_  =  p  j(wt  -  kd') 

Hr  r  » 


where  d'  is  the  reflected  wave  front  vector,  cP  =  -x  cos  a  +  sin  a. 
Morfey  [11]  considered  a  flanged  annular  duct  opening  and  obtained  an 
expression  for  the  impedance  of  the  (m,0)  modes  for  the  case  where 
the  annulus  width  is  small  compared  with  the  wave  length.  This 
expression  for  the  modal  resistance  t  and  reactance  x  is: 

t  =  pv  J  (2mv)  -  ~  [S.p(l  +  v)  -  S.p(l  -  v) ]  +  —  pv2  sin  p 

Zm  7T  1  1  7T 


and 


tt/2 


X  =  pvfi  (2mv)  +  -2~ 
2m  it 


h/4 


cos  2m]>  ,  ,  .  2v  „  .  .  _  , 

- ; - —  dy  +  —  (fn2’sm  >.  +  S.u) 

sin  ij/  tt  v  i  ' 


PV 


2C.P  +  in 1 1  -  v  |  -  C . P ( 1  +  v)  -  C.p(L  -  v) 


where 


and 


kr 


(F3) 


mh 

2r 


P 


with 
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h  =  the  annulus  width, 

2  12  2 

r  =  the  root  mean  square  of  hub  and  tip  radii  r  =  — (r  +r,  ), 

o  o  2  a  b 

J2m  =  integration  of  the  Bessel  function  of  first  kind  of 
order  2m, 

ft2m  =  integration  of  the  Weber  function  of  order  2m, 

=  Sine  integral  function, 

=  Cosine  "integral  function 

and 

in  =  natural  logarithm  function. 

Then,  the  reflection  coefficient  can  be  calculated  by  Equation  (8.29) 
of  Reference  [41] : 

p 

R  =  -i  =  (t  cos  a  -  1)  -  j(x  cos  a)  . 

(x  cos  a  +  1)  -  j(x  cos  a) 


F . 2  Reflection  of  the  Rotor  Blade  Row 
Aimet  [42]  considered  an  obliquely  incident  sound  wave  impinging 
on  a  stationary  blade  row  and  obtained  the  direction  and  amplitude  of 
the  reflected  wave  as: 


R  = 


H  sin  a  (2'j i  sin  0  -  H  sin  a) 

2  2 
(M  sin  a  +  i  ii)  H 


and 


<>  =  tan 


-1  2i[>  sin  0  -  H  sin  a 
2 

2<[i  cos  0  -  H  cos  a 


(F5) 


where  u  and  $  are  the  angle  of  incidence  and  the  angle  of  reflection, 
respectively,  relative  to  the  blade,  and  where 
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M  =  Mach  number, 

4>  =  cos  (0  -  a)  -  M  cos  0, 

2  2  2 

H  =  1  -  M  cos  0 


and 


*•  =  lattice  coefficient. 


Writing  the  incident  wave  as: 

pi  =  p'  ~  ^d) 

i  * 

where  d  =  -x  cos  (ct  -  £)  +  y  sin  (a  -  £)  is  the  wave  front  vector, 
then,  the  reflected  wave  should  be: 

p*  =  RP’.  ei(wC  "  kd,)  , 

r  l  * 

where 


d'  =  x  cos  (>j)  +  £)  +  y  sin  (£  -  £) 


and  5  is  the  blade  stagger  angle.  The  following  summarizes  the 
reflection  considerations  for  two  cases: 

1.  Reflection  at  the  duct  opening:  (a)  change  the  sign  of 
x-component  of  the  wave  number,  (b)  attenuate  the  pressure 
amplitude  by  a  factor  |r|,  and  (c)  shift  the  temporal  phase 
angle  by  the  phase  angle  of  the  complex  quantity  R. 

2.  Reflection  at  the  rotor  blade  row:  (a)  the  angle  of 
reflection  is  not  equal  to  the  angle  of  incidence,  i.e., 
the  helix  angle  changes  after  the  reflection,  (b)  the 
pressure  amplitude  is  attenuated  by  the  factor  R,  and 
(c)  there  is  no  temporal  phase  shift. 

For  the  case  of  the  AFRF  operating  near  the  cutoff  frequency, 
Aimet’s  expression  for  reflection  coefficient  at  the  rotor  blade  row 
(R^)  is  applicable  because  the  wavelength  is  long  (_  4  ft)  compared 
with  the  blade  spacing  (0.45  ft),  and  rotor  chord  length  (0.5  ft),  i.e.. 
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Che  coordinate  is  only  slowly  varying  with  time.  The  restrictions  set 

in  Moriey's  [11]  calculation  for  the  reflection  coefficient  at  the  duct 

opening  (R  ):  kh  <  1,  h  <  r  can  also  be  satisfied  (r  =  0.692  ft, 
o  o  o 

h  =  0.50  ft,  k  =  1.5  ft"1). 

A  computer  program  was  written  to  calculate  the  sound  pressure 
after  R  times  of  reflection  between  r  r.e  two  faces  at  different  measure¬ 
ment  locations  (x,3)  near  the  duct  ill  for  a  known  acoustic  pressure 
source.  Only  the  plane  wave  and  the  first  higher  order  mode  (1,0)  were 
considered.  The  program  reads  in  the  modal  source  pressure  p^  and 
which  can  be  obtained  from  the  program  for  the  unsteady  lift  of  a 
cascade.  The  impedance  of  the  duct  opening  is  a  function  of  frequency, 
and  which  is  calculated  by  Equation  (F4)  varies  with  angle  of  inci¬ 
dence.  In  Equation  (F5) ,  is  a  function  of  the  incident  angle 
relative  to  the  blade  row  for  a  specific  duct  geometry.  The  amplitude 
and  phase  of  the  resultant  pressure  is  then  the  complex  summation  of 
the  incident  and  reflected  waves. 

Plots  of  the  steady  state  pressure  distribution  versus  duct 
axial  length  at  different  operating  conditions  are  shown  in  Figures  FI 
and  F2.  It  can  be  seen  that,  for  the  plane  wave  and  (1,0)  mode 
whose  cutoff  ratio  is  greater  than  1.1,  N  =  1  is  enough  to  reach  a 
steady  state  or  a  convergent  solution.  This  is  because  the  averaged 

reflection  coefficients  R  and  R  are  0.9  and  0.1  respectively,  so 

or 

the  effect  of  the  second  reflection  only  has  a  9  percent  contribution 
to  the  resultant  magnitude  at  most.  However,  for  the  (1,0)  mode  slightly 
above  cutoff,  N  =  3  is  required  to  reach  convergence.  This  occurs 
because  the  second  reflection  at  the  blade  row  is  not  negligible  for  a 


NORMALIZE  PRESSURE, 


NORMALIZED  PRESSURE,  p/p 
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Figure  F2  Steady  State  Wave  Form  After  Reflections  for  (1,0)  Mode 
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helix  angle  near  90  degrees,  i.e.,  the  incident  angle  is  larger  relative 
to  the  rotor  blade  row.  is  a  function  of  stagger  angle,  so  it  varys 

along  the  span.  To  examine  this  s'panwise  variation,  the  convergent 
pressure  distribution  was  plotted  against  stagger  angle  at  hub,  midspan, 
and  tip,  and  it  was  found  that  the  greatest  difference  for  the  cases 
considered  is  8.5  percent.  A  different  spin  direction  relative  to  the 
spin  direction  of  the  rotor  has  a  different  incidence  angle  relative  to 
the  rotor  blade  row.  The  pressure  distribution  for  the  co-spin  (8-cycle 
inlet  distortion)  and  the  counter-spin  (10-cycle  inlet  distortion)  has 
been  compared.  These  two  cases  essentially  result  in  the  same  standing 
wave,  since  R^_  is  very  small  compared  with  R^  so  that,  theoretically, 
spin  direction  does  not  have  a  strong  influence. 


F.3  Modal  Decomposition  with  Reflection 

A  data  reduction  program  was  written  to  calculate  the  source  modal 

pressure  p  and  p,  at  the  rotor  blade  row.  It  has  been  shown  that,  for 
o  1 

all  the  cases  except  the  (±1,0)  mode,  whose  cutoff  ratio  is  slightly 
greater  than  unity,  one  reflection  is  sufficient  to  reach  a  steady  state 
solution  for  the  experimental  conditions.  Thus,  only  the  first  reflec¬ 
tion  is  considered. 

The  resultant  measured  pressure  can  be  written  as: 


(i)  below  cutoff 


m  o 

(ii)  at  cutoff 


-ikx  ,  „  ikx 

p  =  +  R„  P„  e 


o  o 


-ikx  ,  n  -10  .  „  ikx 

p=pe  +P,e  +Rpe 
Km  o  1  o  o 


,  and 
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(iii)  above  cutoff 


-ikx  i(-akx  -  0)  ikx 

p-pe  +P,e  +Rpe 

m  o  1  o  o 


+  R1  P1  C 


i(akx  -  0) 


(F6) 


where  R  and  R.  are  the  complex  reflection  coefficients  at  the  duct 
o  1 

opening  for  the  plane  wave  and  (-1,0)  mode,  respectively. 

The  unknowns  in  case  (i)  are  p  and  R  ,  since  the  (-1,0)  mode  quickly 

o  o 

decays.  For  case  (ii)  ,  there  are  three  unknowns  p  ,  p  1  and  R.,  because 

o  ,  1 

there  is  no  reflection  of  the  (-1,0)  mode  whose  cutoff  ratio  is  smaller 

than  1.0,  so  only  three  measurement  points  are  needed.  For  case  (iii), 

there  are  four  unknowns:  p  ,  p, ,  R  ,  and  R, ,  i.e. ,  at  least  four  measure- 

o  1  o  1 

ment  points  are  required  to  obtain  a  solution.  The  values  of  R^  and  R^ 
obtained  from  the  solution  can  be  used  as  a  check  with  the  theoretically 
predicted  reflection  coefficients. 


APPENDIX  G 


CALIBRATION  OF  UNSTEADY  FORCE  GAGE 

As  discussed  in  Section  2.4.2,  the  frequency  response  of  the  unsteady 
lift  gage  (judging  by  the  broadband  spectrum  analysis)  appeared  to  be 
flat  over  the  frequency  range  of  interest,  with  a  resonance  about  500  Hz. 
The  result  indicated  that  a  static  calibration  of  the  lift  gage  could  be 
employed  over  the  flat  portion  of  the  frequency  response.  To  confirm 
this  flat  response,  a  dynamic  calibration  of  the  blade  segment  was 
attempted  with  a  small  electro-dynamic  shaker  and  a  calibrated  unsteady 
force  gage.  However,  with  these  components  attached  to  the  blade  segment, 
it  was  found  that  the  resonance  frequency  was  lowered  to  about  380  Hz. 

The  shaker  and  force  gage  were  then  detached  from  the  blade  and  the 
voltage  output  was  observed  on  an  oscilloscope  after  tapping  the  blade 
with  the  end  of  a  pencil.  Measuring  the  distance  between  peaks  in  the 
decaying  sine  wave  showed  the  resonant  frequency  to  be  approximately 
500  Hz,  which  indicated  that  the  shaker  assembly  was  responsible  for  the 
difference.  Another  confirmation  that  the  added  mass  lowered  the 
resonant  frequency  of  the  system  was  the  measurement  of  a  resonant 
frequency  of  approximately  450  Hz  when  a  7  gram  accelerometer  was  next 
attached  to  the  blade  and  tapped  as  previously. 

In  order  to  avoid  having  any  added  mass  on  the  instrumented  blade 

segment,  a  seismic  technique  to  measure  the  frequency  response  was 

employed.  The  blade  was  mounted  on  a  large  (50  lb  )  shaker  with  a  7 

rms 

gran  accelerometer  on  the  rigid  portion  of  the  blade,  as  shown  in  Figure 


Gl.  The  input  acceleration  level  was  maintained  constant  for  several 
frequencies  over  the  range  of  interest.  Since  the  acceleration  level 
and  the  mass  of  the  blade  segment  are  constant ,  the  force  on  the  segment 
is  also  constant  by  Newton's  second  law. 

A  voltmeter  was  used  to  set  the  acceleration  level  that  produced 
an  output  voltage  from  the  lift  gage  typical  of  the  values  obtained 
during  testing  in  the  axial  flow  research  fan.  The  frequency  response 
for  two  other  equivalent  force  levels  were  similarly  evaluated  to  check, 
linearity  of  the  lift  gage. 

The  results  of  this  seismic  calibration  are  shown  in  Figure  G2  and 
indicate  a  flat  response  and  linear  behavior  of  the  gage  through  the 
frequency  range  of  interest.  These  results  confirm  that  a  static  cali¬ 
bration  of  the  lift  gage  can  be  employed  and,  in  fact,  for  the  very  low 
mass  blade  segment  employed,  this  is  the  most  accurate  technique  to 
employ.  The  arrangement  used  in  the  static  calibration  is  shown  in 
Figure  G3.  Weights  were  applied  in  the  direction  of  the  lift  force 
by  a  wire  passing  over  a  pulley  mounted  on  3  ball  bearing.  Simultaneously, 
weights  were  added  to  simulate  the  effect  of  the  centrifugal  force  on 
the  instrumented  blade  segment.  Figure  G4  shows  the  results  of  this 


calibration. 


VOLTAGE  OUTPUT  (volt 
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Figure  G3  Instrumentation  of  Static  Calibration  of  Unsteady  Force  Gaye 
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